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A  COMPREHENSIVE  ANALYTICAL  MODEL  OF 
ROTORCRAFT  AERODYNAMICS  AND  DYNAMICS 

Part  III:  Program  Manual 


Wayne  Johnson 

Ames  Research  Center 
and 

Aeromechanics  Laboratory 
AVRADCOM  Research  and  Technology  Laboratories 

SUMMARY 

The  computer  program  for  a  comprehensive  analytical  model 
cf  rotorcraft  aerodynamics  and  dynamics  is  described.  This 
analysis  is  designed  to  calculate  rotor  performance,  loads,  and 
noise;  the  helicopter  vibration  and  gust  response;  the  flight 
dynamics  and  handling  qualities;  and  the  system  aeroelastic  stability. 
The  analysis  is  a  combination  of  structural,  inertial,  and 
aerodynamic  models ,  that  is  applicable  to  a  wide  range  of  problems 
and  a  wide  class  of  vehicles.  The  analysis  is  intended  for  use  in 
the  design,  testing,  and  evaluation  of  rotors  and  rotorcraft,  and 
to  be  a  basis  for  further  development  of  rotary  wing  theories. 

This  report  documents  the  computer  program  that  implements  the 
analysis . 
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The  total  length  of  the  block  is  given  in  parentheses  after  the  label. 
Then  all  variables  in  the  block  are  listed.  The  left-hand  column  gives 
the  variable  name,  and  the  right-hand  column  gives  the  location  cf  the 
variable  in  the  common  block.  Finally,  a  description  of  the  variable  is 
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j  =  1  to  MPSI,  n  =  1  to  max(MHARM,MHARMF*NBLADE) 
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RTR1GM(1070) 


OMEGA 

rotor  speed  n  (rad/sec) 

1 

MTIP 

tip  Mach  ntunher  riR/cg 

2 

GAMMA 

Lock  number 

3 

CMEAN 

mean  chord  c^, 

4 

IB 

characteristic  inertia  I, 

D 

5 

NBM 

number  of  bending  modes 

6 

NTM 

number  of  torsion  modes 

7 

NGM 

zero  if  no  gimbal  or  teeter  mode 

8 

NBHT 

number  of  mean  bending  deflection  modes 

9 

GLAG 

HLD 

®lag  o 

V^b^ 

10 

11 

DZLD 

* 

12 

CGC 

*^GC  ‘^GC'^^^^b^  (or  C^  = 

13 

CGS 

14 

NUGC 

(or  0^) 

15 

NUGS 

'^QS 

16 

CTO 

collective  control  damping  Cq./I^jQ. 

17 

CTC 

cyclic  control  damping 

18 

CTR 

rotating  control  damping 

19 

RA(30) 

aerodynamic  radial  stations,  r^^,  i  =  1  to  MRA 

20 

dra(3o) 

aerodynamic  segment  length  ^r^^,  i  =  1  to  MRA 

50 

FTIP(30) 

tip  loss  multiplicative  factor  f. ,  i  =  1  to  MRA  80 

PSI21H 

(rad),  0.  for  rotor  #1  (for  BODYM, 
M0TNH7  WAKEN,  PNGNM) 

110 

PSI21W 

.  (rad),  — for  rotor  (for 

waken;  wakec) 

11^ 

MUX 

A 

112 

MUY 

113 

MUZ 

114 

RGUST(3,3) 

"^G 

115 

chub(6,i6) 

c 

124 

CBHUB(^3.3) 

~  (including  factor  H  fAl) 

220 

CHUBT(16,6) 

T 

c 

229 
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ALFHP 

PSIHP 

MAT 

CD(2) 

CPSI(2) 

PINTER(36) 

PBURST(36) 


EIXXB(51) 

EIZZB(51) 

MASSB(51) 

TWISTB(5l) 

CENT(5l) 


ITHETB(5l) 

GJB(51) 


MASSI(51) 

ITHETI(51) 

XII(31) 

XCI(51) 

TWISTI(31) 

KP2I(51) 

IPITCH 

KTO 

KTG 

KTR 


(deg) 

'^HP 

\t 

for  drive  train 
for  drive  train  motion 
burst  tip  vortex  in  wake  model 

^inter  J  =  1  to  HPSI 

(rad)  at  j  =  1  to  MPSI 

inertial  and  structural  data  at 
r  =  e  +  (j-l)Ar,  j  =  1  to  HRB+1 

nV/Bi 

o.Vm”' 

'  zz 
m 


inertial  and  structural  data  at 
^  ”  ^FA  j  =  1  to  MRB+1 

xiV/gj 


inertial  data  at 
r  =  (j-l)Ar,  j  = 

mnVi^ 

Xj/R 

x^/H 


k  J/ 


1 


to  MRM+1 


blade  pitch  inertia  (slug-ft^  or  kg-m^) 


control  system  stiffness 
(ft-lb/rad  or  m-N/rad) 
collective 
cyclic 

reactionless 


RTRICH 

325 

326 
32? 
328 
330 

332 

368 


404 

455 

506 

557 

608 


659 

710 

761 

812 

863 

914 

965 

1016 

106? 

1068 

1069 

1070 
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RH1CH( 12792) 
HRTR( 16, 16,21) 

HB0DY(l6,6,10) 

HENG(6,10) 


complex  rotor  transfer  function  matrix,  ; 
size  NBM+NTM+NGM;  n  =  0  to  MHARM  ” 

complex  airframe  transfer  function  matrix, 

;  n  =  p  =  1  to 

MHARMF 

complex  drive  train  transfer  function  matrix, 

"n  ‘  <=D  P  -  1  to 

MHARMF 


1 

10753 

12673 


-26- 


B0DifCM(446) 

AMO DEI (6, 16) 

rotor  HI  hub  (dimensionless) 

1 

AH0DS2(6,16) 

(  '^k*  ^jj)  at  rotor  ifZ  hub  (dimensionless) 

61 

KHSTCl(lO) 

pitch/mast-bending  coupling  (dimensionless) 

Kj^Cj^  rotor  #1 

121 

KHSTSl(lO) 

K^Sj^  for  rotor  #1 

131 

KMSTC2(10) 

K|^Ck  rotor  liZ 

141 

KMSTS2(10) 

Kmsjj  for  rotor  #2 

151 

ADAMPA(lO) 

aerodynamic  damping  (2X/’r*aA)(q/v)Fn  ^ 

'ik^k 

l6l 

ACNm(4,10) 

control  derivatives  (2V/v-aA)  q  F^^^ 

171 

AHASS(IO) 

* 

\ 

211 

ADAMPS(lO) 

"kSs^^k 

221 

ASPRNG(lO) 

<  -k 

231 

MSTAR 

m" 

241 

MSTARG 

M^'g 

242 

istar(3,3) 

I* 

243 

CWS 

C„/^  =  (a/2  2()  H*g 

252 

HHASS 

aircraft  mass  (slug  or  kg) 

253 

NAM 

number  of  airframe  modes 

254 

RSF10(3,3) 

aircraft  description  (  ^^  =  =  O) 

Rgp  for  rotor  #1 

255 

RSF20(3,3) 

Rgp.  for  rotor  ^Z 

264 

RHUB10(3) 

"x  at  rotor  UX  hub 

273 

nKU320(3) 

r  at  rotor  #2  hub 

276 

RWBO(3) 

'r  at  wing/body 

279 

RHTO(3) 

r  at  horizontal  tail 

282 

RVTO(3) 

r^  at  vertical  tail 

285 

ROFFO(3) 

r  off  rotor 

288 

BODYCM 


aircraft  description 


HSF1(3,3) 

Rgp  for  rotor 

291 

RSF2(3,3) 

Rgp  for  rotor 

300 

RHUBl(3) 

"t  at  rotor  li\  hub 

309 

RHUB2(3) 

Ir  at  rotor  tt2.  hub 

312 

RWB(3) 

at  wing/body 

315 

RHT(3) 

r  at  horizontal  tail 

318 

RVT(3) 

r  at  vertical  tail 

321 

R0FF(3) 

^  off  rotor 

324 

TCFE(11,5) 

"^GFE 

327 

VXRSKF(3) 

(^x  )  R^  kp 

382 

MVXRE(3,3) 

-M*  (Vx  )  R^ 

385 

GMTRX(3,3) 

G 

394 

IBODY(3,3) 

R*^  I*  R 
e  e 

403 

REULER(3,3) 

R 

e 

412 

RFV(3,3) 

^FV 

421 

RFE(3,3) 

^FE 

430 

KE(3) 

439 

VELF(3) 

V 

^(42 

VCLIHB 

^climb 

445 

'side 

446 

ENGNCH(131) 


QTHRTL 

^  * 

1 

lENG 

2  T  * 

^E  ^E 

2 

KMIl 

3 

KMI2 

%2 

4 

KHR 

^HR 

5 

KMEl 

6 

KME2 

^ME2 

governor  proportional  gains, 

7 

KPGOVE 

engine 

8 

KPGOVl 

rotor  #1 

9 

KPGOV2 

rotor 

10 

NDM 

number  of  drive  train  modes 

governor  time  lag, 
engine 

11 

TIGOVE 

12 

TlGOVl 

rotor  #1 

13 

T1GOV2 

rotor  #2 

governor  time  lag,  x 

14 

T2GOVE 

engine 

15 

T2GOV1 

rotor  #1 

16 

T2GOV2 

rotor  #Z 

17 

QEDAMP 

9  * 

18 

IRSTAR 

I* 

^  -1 
mass  matrix  for  H 

n 

19 

HENG(6,6) 

20 

SENG(6,6) 

-1 

spring  matrix  for  H 

56 

DENG(6,6) 

-i 

damping  matrix  for 

92 

HENG0(2,2) 

Hq  ^  for  static  elastic  motion 

128 

GUSTCM( 12989) 

VGWBV(3) 

gust  components,  velocity  axes 
at  wing/body, 

1 

VGHTV(3) 

at  horizontal  tail,  gj^ 

4 

VGVTV(3) 

at  vertical  tail,  gy 

7 

VGR1V(3,30,36) 

at  rotor  ^1,  gir^, 

10 

vgr2V(3,30,36) 

at  rotor  ^2,  gir^, 

3250 

VGHUB1(3) 

at  rotor  hub,  "g  (for  wake  geometry) 

6490 

VGHUB2(3) 

at  rotor  ffZ  hub,  "g  (for  wake  geometry) 

6493 

VGWBF(3) 

gust  components,  F  axes 
at  wing/body,  "g^ 

6496 

VGHTF(3) 

at  horizontal  tail,  gj^ 

6499 

VGVTF(3) 

at  vertical  tail,  ^ 

6502 

VGR1S(3,30,36) 

gust  components,  S  axes 

at  rotor  #1, 

6505 

VGR2S(3,30,36) 

at  rotor  #2,  g(rj^, 

9745 

VPTRAN(5) 

transient  control 

A'vp  =  (  Sq  ■^s  ^  ^t^*^ 

12985 
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C0NTCM(32) 

VCNm(ll) 

control  vector  (rad):  „ 

1 

Totorffl  TotoiffZ  airframe 

THETFT 

12 

PHIFT 

^FT 

13 

THETFP 

©^p  (rad) 

lit 

PSIFP 

Vpp  (rad) 

15 

THETAT 

(rad) 

16 

PSIT 

(rad) 

17 

DVB0DY(6) 

airframe  motion  (dimensionless) 

( ‘ft'  ft  'Vp  ^  ft  Zp) 

18 

DOHEGA 

'Vg  (static;  dimensionless) 

24 

DDZF 

*ft  (dimensionless) 

25 

vpilot(5) 

pilot  control  vector  (rad): 

Vp-(i,  ^3  ip  i/ 

26 

TGOVRl 

TGOVR2 

^  ^  ^go  vr  ^  ro  tor^^l  ^  ^ 

^'^‘^govr^rotor#2 

31 

32 

-31- 


CONVCH(80) 

mean  square  motion  (rotor  n^l) 

BIHS(IO) 

1 

T1MS(5) 

€) 

11 

BGIMS 

16 

PlHS(l6) 

4> 

17 

PS1HS(6) 

mean  square  motion  (rotor  ^2) 

33 

B2HS(10) 

39 

T2MS(5) 

& 

'49 

BG2MS 

54 

P2MS(l6) 

55 

PS2MS(6) 

t 

71 

GIHS 

mean  square  circulation  (rotor  n^l) 

77 

G2MS 

mean  square  circulation  (rotor  #2) 

78 

COUNTM 

integer  parameter:  ntunber  of  motion  iterations 

79 

COUNTC 

integer  larameter:  number  of  circulation 
iterations 

80 

-32- 


T750LD 

old  (initialized  to  1000.) 

1 

NBMCLD 

old  NBM  (initialized  to  O) 

2 

NTOOLD 

old  NTM  (initialized  to  O) 

3 

mj(20) 

bending  frequency  ,  i  =  1  to  NCOLB  (per  rev) 

4 

NUNR(20) 

nonrotating  bending  frequency 
i  =  1  to  NCOLB  (rad/sec) 

bending  mode  displacement 
at  radial  station  r  = 

0 

^NRi* 

i  =  1  to  NBM, 

24 

ETA(2,10) 

^FA 

44 

ETA(2,10) 

^PB 

64 

ETA(2,10) 

^ROOT 

84 

ETA (2, 10) 

1 

104 

ETA(2,10.11) 

(j-l)O.l,  j  -  1  to  11 

124 

ETA(2,10,51) 

(.1-1)^.  j  =  1  to  MRM+1 

344 

ETA(2,10,30) 

r 0=1  to  HRA 

bending  mode  s].ope  'n i  =  1 
station  r  =  ' 

to  NBM,  at  radial 

1364 

ETAP(2,10) 

^FA 

1964 

ETAP(2,10) 

’^PB 

1984 

ETAP(2,10) 

^ROOT 

2004 

ETAP(2,10) 

1 

2024 

ETAP(2,10,11) 

(o-l)O.l,  j  =  1  to  11 

2044 

ETAP(2,10,5l) 

(j-l)Ar,  0  1  to  MRM+1 

2264 

ETAP(2:10,30) 

r . ,  0  =  1  to  MRA 

J 

bending  mode  curvature  i 

at  radial  station  r  = 

=  1  to  NBM, 

3284 

ETAPP(2,U) 

^FA 

3884 

ETAPP(2,iO) 

^PB 

3904 

ETAPP(2,10) 

^ROOT 

3924 

ETAPP(2,10) 

1 

3944 

ETAPP(2,10,ll) 

(o-i)O.l,  j  =  1  to  11 

3964 

ETAPP(2,10,5l) 

(’-l)Ai^i  0  =  1  to  MRM+1 

4184 

ETAPP(2,10,30) 

r . ,  j  =  1  to  MRA 
<3 

5204 

-33- 


ETAPH(2,10) 

WT{11) 


ZETA(5.11) 

ZETA(5.51) 

2ETA(5,30) 


zbtap(5,u) 

ZETAP(5,51) 

ZF,TAP(5,30) 

KPB(IO) 


bending  mode  slope  at  hinge, 

1  -  1  ti  KCOLT+l, 

control  system  frequency  (per  rev) 
collective 
cyclic 

reactionless 

torsion  mode  displacement  i  =  i  to  NTW 
at  radial  station  r  =  i  to  nto, 

(j-l)O.l.  j  =  1  to  11 
(j-l)Ar,  j  =  1  to  MRM+l 
Ty  j  =  1  to  HRA 

iSlon  7t"  1  “  1  to  »™,  at  radial 

(j-i)O.l,  j  =  1  to  11 
(j-l)Ar,  i  =  1  to  MRM+l 
rj,  j  =  1  to  MRA 

pitch/bending  coupling  Kp^,  i  =  1  to  NBM 
pitch/gimbal  coupling  Kp^ 

^^rA,  (fad) 

(rad) 

(rad) 

ScA_  (rad) 


HDICM 


6298 

6353 

6608 

6738 

6768 

6769 

6770 

6771 

6772 

6773 


CTOLD 

old 

1 

1 

CMXOLD 

old 

2 

CMYOLD 

old 

3 

gamold(30,36) 

old  (i  =  1  to  MRA,  j  =  1  to  MPSI) 

4 

CRC0LD(36) 

old  max  (j  =  1  to  MPSi) 

1084 

vind(3,30,36) 

V(r,  in'.)  (i  =  1  to  MRA,  j  =  1  to  MPSI) 

^  J 

1120 

LAMBDA 

'^tpp 

4360 

FGE 

^GE  ""  ^  (cos£/4z)^  (l.  if  OGE) 

4361 

COSE 

cosC 

4362 

7AGL 

^AGL 

4363 

vint(3,30,36) 

VintC^i.^j)  (i  =  1  to  MRA,  j  =  1  to  MPSi) 
at  other  rotor 

4364 

vorh(3,36) 

\i^t('^j)  (o  =  1  to  MPSI),  at  other  rotor  huh 

7604 

LAMBDI 

mean  at  other  rotor 

-A  int 

7712 

vwb(3,36) 

(j  *=  1  to  MFSI),  at  wing/body 

7713 

vht(3,36) 

(o  1  to  MPSi),  at  horizontal  tail 

7821 

vvt(3,36) 

(j  =  1  to  MPSI),  at  vertical  tail 

7929 

voff(3,36) 

'\Q(Vj)  (j  =  1  to  MPSI),  off  rotor  disk 

3037 

LAMBDW(3) 

mean  \  at  wing/hody 

3145 

LAMBDH(3) 

mean'^u,  at  horizontal  tail 

3148 

LAMBDV(3) 

mean  \y,  at  vertical  tail 

3151 

LAMBDO(3) 

mean  off  rotor  disk 

3154 

T?TMrniI^ 

UXil  Xlf  \  ^  / 

=  V  G  Vr:ip)/('CSR') 

'll  'ref 

3157 

einth(3) 

31 60 

EINTV(3) 

"  K,0„BgI  (-•«g)(IXR)/(SXR)^3f 

8163 

Inertia  coefficients 


INC1CM(4365) 


MB 

SB 

10 

IQ(IO)  ^ 

SQ(2,10) 

IQA(2,10) 

IQDQ(10,10) 

IQDQT(10,10,4) 

IQDP(IO) 

IQDPT(10,4) 

IQDB(IO) 

IQDBT(10,4) 

SQDDP(10,5) 

SQDDPT(10,5,4) 

SQP(10,5) 

SQPT(10,5,4) 

IQO(IO) 

IQ0D^(2,1C) 

IQ0DQT(2,10,4) 

IQODP 

IQ0DPT(4) 

IQODB 

IQ0nBT(4) 

SQ0DDP(2,5) 

SQ0DDT(2,5,4) 

IFXO 

IMXO 

IP(5) 

IPA(2,5) 

IPAT(2,5.4) 

SP(2,5) 

SPT(2,5,4) 

ipddp(5 ,5) 

IPDDPT(5'i5A) 

ipddtt(5,5A.^) 

IFP(5,5)  ^ 

SPDDQ(5.10) 

SPDDQT(5,10,4) 

IP0(5) 

spq(5.10) 

SPQT(5,10,4) 

XAPQ(2,5A.30) 

mqdq(io,io) 

MQr©(lO) 

MQP(10,5) 

mdq(io) 

MDB 

MP(5) 


Xj^j  at  i  =  1  to  MBA 
Aerodynamic  spring  and  damping 


1 

2 

3 

4 
14 
34 
54 

154 

554 

564 

604 

614 

654 

704 

904 

954 

115^^ 

1164 

1184 

1264 

1265 

1269 

1270 

1274 

1284 

1324 

1325 

1326 
1331 
1341 
1381 
1391 
1431 

1456 

1556 

1956 

1981 

2031 

2231 

2236 

2286 

2486 

3686 

3786 

3796 

3846 

3856 

3857 
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INCICM 


QDZ  3B62 
QT  3863 
MPIXi(5,10)  3864 
MPDB(5)  3914 
MPDP(5.5)  3919 
MPP(5.5)  3944 
IQDQS(l0,10)  Inertia  coefficients,  summed  over  q.  3969 
IQDPS(10}  ^  4069 
IQDBsho)  4079 
SQDDPS(10,5)  4089 
SQPS(10,5)  4139 
IQ0D^S(2,10)  4189 
IQODPS  4209 
IQODBS  4210 
SQ0DDS(2,5)  4211 
IPAS(2,5)  4221 
SPS(2,5)  4231 
IPDDPS(5,5)  ^241 
SPDDQS(5,10)  4266 
siQS(5,io)  4316 


(NBH=10,  NTM=5,  NBM'1'=4,  mra=3o) 
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WKV1GM(’816‘5^ 

CTOLD 

old 

1 

CMXOLD 

old 

2 

CMYOLD 

old 

3 

GAM0LD(30,36) 

ol'i  (i  "  1  to  MRA,  j  =  1  to  MPSI) 

4 

CRCOLD(36) 

^Id  max  (j  =  1  to  MPSi) 

1084 

vind(3,30,36) 

(i  =  1  to  MRA,  j  =  1  to  MPSI) 

1120 

LAMBDA 

mean  V , 

tpp 

4360 

FGE 

^GE  "  "1-  (cos£/4z)2  (i.  if  qgE) 

4361 

COSE 

COSC 

4362 

ZAGL 

fAGL 

4363 

vint(3,30,36) 

'^int^’^i*  (i  =  1  to  MRA,  j  =  1  to  MPSi) 

4364 

V0RH(3,36) 

at  other  rotor 

'^int^'^j^  (j  =  1  to  MPSi),  at  other  rotor  hub 

7604 

LAMBDI 

jnean  at  other  rotor 

7712 

vwb(3,36) 

(j  =  1  to  MPSi),  at  wing/body 

7713 

vht(3,36) 

(j  =  1  to  MPSi),  at  horizontal  tail 

7821 

vvt(3,36) 

=  1  to  MPSi),  at  vertical  tail 

7929 

voff(3,36) 

(i  “  1  to  MPSi),  off  rotor  disk 

3037 

LAMBDW(3) 

mean  \y,  at  wing/body 

3145 

LAMBDH(3) 

mean’^VTi.  at  horizontal  tail 

3148 

LAMBDV(3) 

mean  at  vertical  tail 

3151 

LAMBD0(3) 

mean  off  rotor  disk 

3154 

EINTW(3) 

%  -  w|.  (-tg)(nR)/(£«)^^j, 

®H  "  Vh^SF 

3157 

EINTH(3) 

8160 

EINTy(3) 

81 63 
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HNH1CM(462) 

ALF(10,6) 

complex  (p  =  1  to  MHARMF),  without  Euler 

1 

angle  contributions 

DALF(10,6) 

complex  (p  ~  1  to  MHARMF) 

121 

DDALF(10,6) 

complex  o<pj^  (p  =  1  to  MHARMF) 

241 

PSIS(IO) 

complex  (P  “  1  to  MHARMF) 

361 

TGOVR(lO) 

complex  ^  MHARMF) 

381 

TMAST(2l) 

complex  (n  -  1  to  MHARM) 

401 

ALF0(6) 

*^static 

443 

DDAL0(6) 

• 

*^static 

449 

DDALFO(6) 

^static 

455 

PSISO 

static 

461 

DPSISO 

static 

462 
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1 


state(30,36,3) 


dclhax(30,36) 

r)CDMAX(30,36) 

dcmmax(3o,36) 

meff(30,36,3) 

AEFF(30,36,3) 

dclds(30,36) 

dgdds(3o,36) 

dgmds(30,36) 

SAVE(30,36,19) 


integer  parameter  defining  stall  state  for 
lift,  drag,  moment  (initialized  to  zero) 


peak  dynamic  stall  vortex  loads  (initialized 
to  zero) 


^^max 

3241 

^^max 

4321 

^"'max 

5401 

effective  environment  for  lift,  drag,  moment 

Mach  number 

6481 

angle  of  attack  oc’g^^ 

9721 

dynamic  stall  vortex  load 

12961 

14C41 

Ws 

15121 

section  aerodynamic  data 

J  K  i>  W 

16201 

(1) 

Up 

(11) 

% 

(2) 

(12) 

°d 

(3) 

(13) 

°m 

(4) 

U 

(14) 

^dradial 

(5) 

e 

(deg) 

(15) 

V^'^m 

(6) 

(deg) 

(16) 

F  Vac 
r  m 

(7) 

CK 

(deg) 

(17) 

z  m 

(8) 

H 

(18) 

M  /ac 

in 

(9) 

cosA. 

riQ^ 

F  /ac 

T'  m 

(10) 

<kc/V 

aerodyneunic  data  at  (r, ,  H^.)  on  disk, 
i=ltoMRA,  j=l  to^HPS^ 


MNR1CM(1112) 

(0 

complex  (i  =  1  to  NBM,  n  =  0  to  MHARH) 

BETA(21,10) 

1 

THETA(21,5) 

complex  (i  =  1  to  NTM,  n  =  0  to  MHARM) 

421 

BETAG(21) 

complex  ^  (n  =  0  to  MHARM) 

631 

phi(io,i6) 

J''*  / 

complex  (i  =  1  to  NAM,  p  =  1  to  MHARMP) 

673 

PSID(10,6) 

co.pi«x  ( A6t  AAgj  -aeg2)pN 
(p  =  1  to  MHARMP) 

993 

-40- 


HNSCM( 12; 

QSSTAT(IO) 

(isk'sutic  C'  “  7 

elastic 

PISTAT 

static 

elastic 

PESTAT 

^^e^static 

elastic 


-41- 


AEF1CM(1548) 

FORCE( 16,36) 

'^Plast  rev.  J 
(dimension  NSM+NTO+NGM) 

FHUB(6,36) 

hub  reactions  (without  rotor  mass  terms) 
F  =  (  ^2Cjj/*<»-a,  iS2C^/^a.,  ^2C,p/-<-a, 

^2CM^’<3-a,  1SZ0l^^^^a.,  -‘XZC^ma.) 

T0RQUE(36) 

SAVE(36,20) 

integrated  aerodynamic  forces 

(i.)-(lO)  Mq^g^pj.Q/ac 

(11)-(15)  Mp^ag^o/ac 

(16)  V/«5-a 

f 


I 


QR1CN(1139) 


QRTR(6) 

rotor  generalized  force,  "q  =  c'^F 

1 

FHUBH(6) 

mean  hub  reaction 

7 

F  =  (^2Cjj/r-a,  X2C,pAr~a, 

1 

^2CMyv-a,  ^S2Ci^y'<7-a, 

for  trim 

CLS 

Cj^/t"  (wind  axes) 

13 

CXS 

(wind  axes) 

14 

CTS 

V- 

15 

CYS 

Cy/g- 

16 

CPS 

Cj/<r- 

17 

for  inflow 

CT 

Ct 

18 

CHX 

CMx 

19 

CHY 

C,. 

20 

4>Iy 

for  trim 

BETAO 

Po 

21 

BETAC 

(ic 

22 

BETAS 

Ps 

23 

for  inflow 

circulation  (i  =  1  to  MRA,  j  =  1  to  MIJl) 

maximum  circulation  P  (j  =  1  to  HPSi) 

<3 


gah(30,36) 

circ(36) 


24 

1104 


OS  per  -9) 

qv;b(6) 

QHT(6) 

qyi'iS) 

SAVE(3l) 


wing-body  generalised  forces 

horizontal  tail  generalized  forces 

vertical  tail  generalized  forces 

airframe  aerodyneimic  data  _ 

(l)  (D/q),,„  ft'^  or  m ' 


(1) 

(D/q),B 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(D/q)^^, 

(«) 

(L/q)HT 

(^) 

(0/q)yY 

(10) 

(l*/  9)y<p 

3  3 

•f*  4-'^  ir\“^ 


„.2  2 
ft  or  m 


(11) 


(1^0 

(15)  € 

(16)  0- 
(17-1^)  t,3 
(20-22) 
(23-25) 

(26-28)  w 

(79) 

(30)  q^^ 

(31)  q„,p 


'sec  or  m/sec 


rad/sec 

dimensionless 


WG1CH(7948) 

RBR(3,36) 

RBT(3,36) 

1 

lOP 

MUTFP(3) 

A'tpp 

21? 

DZT(l44) 

prescribed  wake,  tip  vortices 

D^(k),  k  =  1  to  KRWG 

220 

DnT(ivO 

D  (k),  k  =  1  to  KRWG 

JfA 

K2T 

K. 

4. 

508 

DZSld^O 

prescribed  wake,  sheet  inside  edge 

509 

DRSl(i44) 

n  (k),  k  =  1  to  KRWG 

653 

K2SI 

^2 

?Q7 

prescribed  wake,  sheet  outside  edge 

D^(k),  k  -  1  to  KRWG 

798 

DRSC(l44) 

D  (k),  k  -  1  to  KRWG 

942 

K2SO 

K^ 

1086 

DFWG(3,2304) 

free  w^e,  tip  vortices 

D(n),  n  =  1  to  KRWG^nPSI 

1087 

n  =  (51-1  )KFV;G  +  k 

((k  =  1  to  KFWG),  ^  =  1  to  HPSI) 
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WKC1CM( 120007) 

t-!R 

total  number  of  points  in  flow  field  at  v;hich 
nonuniform  induced  velocity  calculated  for 
each  azimuth  (ML+HI+MVJ+MJ{i-MV+r'iO) 

1 

ML 

number  of  points  on  this  rotor 

(MRL  if  INFLOW(l)  =  I;  zero  otherwise) 

2 

MI 

number  of  points  on  other  rotor 
(MRL  of  other  rotor  if  INFLOW (2)  =  3; 

I  if  INFLCW(2)  =  2;  zero  otherwise) 

3 

MW 

ntunber  of  points  on  wing-body 

(l  if  INP’L0W(3)  =  2;  zero  otherwise) 

4 

MH 

number  of  points  on  horizontal  tail 
(l  if  INFLCW(4)  =  2;  zero  otherwise) 

5 

MV 

number  of  points  on  vertical  tail 
(l  if  IN?LCW(5)  =  2;  zero  otherwise) 

6 

MO 

number  of  points  off  rotor  disk 
(l  if  INFLCW(6)  =  I;  zero  otherwise J 

7 

C( 3, 20000) 

"0(0),  n  =  I  to  HPSI^^MR*MPSI 

8 

CN'W(3,20000) 

4'/n^),  nj^  =  1  to  N?,G^(KmHl)W*MPSI 

60008 

"NwWlW^ 

n  =  ((S^  -  l)''MR  +  k  -  l)*MPSI  +  j 

r  -i  =  1  +0  KWT^  V  =  1  +r.  9  =  1  -rr.  MP^T  ^ 


'•••/  »  <\ 


"fT//  "  ■*■  ^"  l)*(KM+l)  +  j-Sl  +KN^rf)*MRG  +  i 

((((i  =  1  to  MRG),  j  =  -il-  KNW  to  ) , 
k  =  I  to  I-IRL),  A  =  1  to  MPSI) 
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0 


AEMNGH(78) 


Q(10) 

qj^ ,  k  ^  1  to  NBM 

1 

DQ(lO) 

n 

DKi(lO) 

21 

P(3) 

Pj^,  k  =  1  to  NTO  (pq  =  P^+  P^) 

31 

DP(5) 

36 

DDP(5) 

41 

PD 

Pd 

46 

DPD 

h 

47 

DDPD 

Pd 

48 

PR 

P^ 

49 

DPR 

Pr 

50 

DDPR 

'4 

51 

BG 

f^G 

52 

DBG 

(^G 

53 

DDBG 

^G 

54 

AHUB(6) 

(without  Euler  angle  contributions  to  oi^) 

55 

DAHDB(6) 

=  (x.  Yv,  Z.  cK  ci  ) 

^  h  ''  h  h  X  y  z ' 

61 

ddahub(6) 

67 

PS 

s 

73 

DPS 

s 

74 

DDPS 

•^s 

75 

fVM 

j,  1  • 

mast- bend 

76 

TG 

^®SOvr 

77 

DTT 

®G  ^  ^G  ^^G 

78 
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LDMNCM(2932) 

SAVEM( 36, 78) 

motion  at  H>j,  j  =  l  to  MPSI 

(refer  to  common  block  AEMNCH  for  contents) 

1 

MB 

SB 

inertial  coefficients  for  section  loads 

2809 

10 

2810 

SQ(2,10) 

2811 

IQA(2,10) 

2812 

IQDQ(2,10) 

2832 

IQDB 

2852 

IQDP(2) 

2872 

SQDDP(2,5) 

2873 

SQP(2,5) 

2875 

IFXO 

2885 

IMX0(2) 

2895 

IPDDP(5) 

2896 

IPP(5) 

2898 

IPA(2) 

2903 

SPDDQ(IO) 

29O8 

SKl(lO) 

2910 

SP(2) 

2920 

IPO 

2930 

2932 

-48- 


FLMCM( 21928) 


A2(6400) 

A1(6400) 

1 

6401 

A0(6400) 

12801 

B(2320) 

19201 

D0F1(80) 

21521 

NAMEX(80) 

21601 

NAMEy(29) 

21681 

MX 

21710 

MXl 

21711 

MV 

21712 

••!G 

21713 

DC'riS(46) 

symmetric  matrices 

21714 

NAMEXS(46) 

21760 

namevs(i6) 

21806 

MXS 

21822 

MX  IS 

21823 

M''S 

21824 

MGS 

21825 

rX)FlA(43) 

antisymmetric  matrices 

21826 

IT  4  4  /  1.^  \ 

i'JrtnSiArt  ) 

21869 

NAHSVA(13) 

21912 

MXA 

21925 

MXl  A 

21926 

M'M 

21927 

MGA 

21928 

variables  (80) 

^  Xo  Mi  Ad+  ) 

K1  KZ  b  e  “  t  govri  *^^govTp' 

controls  (29) 

''  °  ("ni  'rz  ’'s  ''p 


A2(30,30) 

^‘2 

1 

Al(30,30) 

901 

A0(30,30) 

*0 

1801 

AA2(30,6) 

2701 

AA1(30,6) 

2881 

aao(30,6) 

3061 

B(30,8) 

B 

3241 

BG(30.3) 

3481 

G2(6,30) 

c. 

3571 

c1(6,30) 

3751 

C0(6.30) 

Cq 

3931 

ca2(6,6) 

^2 

4111 

CA1(6,6) 

4147 

GA0(6,6) 

4183 

DG(6,3) 

\ 

4219 

variables  (30): 
controls  ( 8) : 
gust(3);  g 
hub  motion  (6); 
hub  forces  ( 6) :  F 


FLMACH(912) 


A2(i6,16)  1 

A1(i6,16)  257 

A0(l6,l6)  aQ  513 

B(l6,4)  b  769 

BG(l6,3)  b^  833 

BL(l6,2)  b^  881 


variables  (I6):  Xg 
controls  (4) :  Vg 
gust  (3)s  g 
inflou(2);  (  ) 
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FLINCMr47?^ 

MASSB 

10 

1 

IQ(IO) 

2 

SQ(10,2) 

3 

IQA(10,2) 

13 

IQDQdO.lO) 

33 

IQDP(IO) 

53 

IQDB(IO) 

153 

SQDDP(10,5) 

163 

SQP(10,5) 

173 

IQ0DQ(10,2) 

223 

SQ0DDP(5,2) 

273 

IP(5) 

293 

IPA(5,2) 

303 

SP(5.2) 

308 

IPDDP(5,5) 

318 

IPP(5.5) 

328 

SPDDQ(5,10) 

353 

SPQ(5,10) 

378 

428 

-52- 


V 


r 


FLAECH(646) 

MQU(lO) 

MQDZ,(in'> 

MQZ(lO) 

MQL(IO) 

MQDB(IO) 

MQB(lO) 

MQDQ(t0,10) 

MQQdO.lO) 

MQP(10,5) 

MMU 

HDZ 

MZ 

ML 

MHB 

MB 

MDQ(IO) 

MQ(IO) 

MP(5) 

TU 

TDZ 

TZ 

TL 

TDB 

TB 

TDQ(IO) 

TQ(lO) 

TP(5) 

HU 

HDZ 

HZ 

HL 

HDB 

HB 

HDQ(lO) 

HQ(10) 

HP(5) 

QU 

QDZ 

QZ 

QL 

QDB 

QB 

QDQ(IO) 

QQ(IO) 

QP(5) 

HR 

RU 

RDZ 

RZ 


-53- 


FLAECM 


RL 

RDB 

RB 

RDQ(IO) 

RQ(IO) 

RP(5) 

MPU(5) 

MPDZ(5) 

MPZ(5) 

HPL(5) 

MPDB(5) 

MPB(5) 

MPEQb.lO) 

MPQ(5.10) 

MPP(5,5) 

MPDP(5,5) 


439 

i440 

441 

442 

452 

462 

46? 

472 

477 

482 

487 

492 

497 

547 

597 

622 


-54- 


(both  rotors  for  flutter  case) 


STDGM(882) 

DERIV(7,21) 

’W'R1(7,21) 

DHVR2(7,21) 

DR''WB(7,21) 

DR'‘HT(7,21) 

DRVVT(7.21) 


1 

148 

295 

442 

589 

736 


variables  (2l): 

( Zp  Xp  Yp  Zp, 

®0  «lc  1,  ®ls  4  i 


equations  (7)s 

(L  M  N  X  Y  Z  Q) 
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S'IWCH(340) 

A2FD(7,7) 

A1FD(7,7) 

A0FD(7,7) 

BFD(7,19) 

TX)FFD(7) 

CCNFD(l6) 

GUSFD(3) 

D0F1FD(7) 

NAMXFD(7) 

NAMVFD(19) 

MXFD 


1 

50 

99 

148 

281 

288 

304 

307 

31i^ 

321 

340 


variables  (7)s 

( F  %  %  ^  ^F  ’^F  *^8^ 

controls  (19)5 

(  a  &  0  ^  ©  © 

'‘^0  Ic  Is  0  Ic  is 

^a  "^r  ^0 

\  '^G  ”g^ 

gust  components  in  wind  axes 


-56- 


TRANCM(62) 

qtrim(6) 

trim  generalized  force  (total) 

1 

CQSTl 

trim  —^ZC^/^a.  (rotor  ffi) 

7 

CQST2 

trim  -32C^5r-a  (rotor  #2) 

8 

IB0DYI(7.7) 

inverse  of  body  inertia 

9 

DCSAS 

SAS  i 

c 

5fl 

DSSAS 

SAS 

59 

TTGOV 

transient  governor 

60 

TIGCA’ 

transient  governor  (^&  )  .  ^ 

govr  rotor//l 

61 

T2GOV 

transient  governor  (aA  )  , 

^  govr '  TO  tOTffl 

62 

-57- 


2.  SUBPROGRAM  FUNCTION  AND  COMMUNICATION 


This  section  describes  the  functions  of  the  subprograms  that 
constitute  the  computer  program.  The  communication  of  the  subprograms 
with  each  other  is  also  described,  in  terms  of  the  input  and  output 
variables,  '^ns  description  begins  with  the  subprogram  name,  and  its 
arguments.  Next  there  is  a  statement  of  the  principal  function  of  the 
subprogram,  and  usually  a  general  reference  to  a  section  in  the  analysis 
development.  Then  notes  about  the  program  content  are  given,  including 
references  to  sections  in  the  analysis  development  as  appropriate. 
Finally  all  the  input  and  output  variables  of  the  subprogram  are  listed . 
The  left-hand  column  gives  the  variable  name  in  the  subprogram,  and  the 
right-hand  column  gives  the  label  of  the  common  blcck  in  which  the 
variable  is  located.  Some  description  of  the  variable  may  be  given  as 
well.  Only  the  subprograms  for  rotor  are  described;  the  subprograms 
for  rotor  -iZ  have  identical  functions  and  structure. 
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MAIN 


Name:  MAIN 

Function:  primary 

job  and  analysis  control 

General  reference: 

section  5»3*5 

CPRTR2 

IDENT(4) 

TRIMCM 

ANTYPE(3) 

FTLEID(4) 

TMDATA 

RESTRT 

JCASE 

TASK 

JOB 

RSWRT 

NCASES 

BLKDAT 

RDFILE 

START 

CASECM 

-59- 


TIMER 


Na..ie:  TIHER(N,I,T) 

Function:  program  timer 

N  integer  parameter  controlling  timing  calculations 

0  initialize 

1  start  timer 

2  stop  timer 

3  print  times 

other  return  present  time 

I  timer  number 

1  case 

2  TRIM 

3  FLUT 

4  STAB 

5  TRAM 

6  STABL 

7  FLUTL 

R  WAKECl ,WAKEC2 
9  GE0MRl,GE0hR2 

10  RAHF 

11  M0DE1,M0DB2 

12  M0TNR1,M0TNR2 

13  PERF 

14  LOAD 

T  elapsed  CPU  time  (sec) 

DEBUG  integer  parameter:  print  time  T  if  GE  1  TMDATA 

ITIB 

IDB(23) 
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INPTN 


Name :  INPTN 


Function: 

input  for  new 

job 

JCA3E 

DLKDAT 

RDFILE 

CASECM 

DEBUGI 

OPREAT)(lO) 

NROTCR 

integer 

parameter:  debug  print  control 

TMDATA 

IXX 

BDDATA 

lYY 

IZZ 

IXY 

IXZ 

lYZ 

ATILT 

FSGG 

BLCG 

WLCG 

WEIGHT 

FILEID(4)  TMDATA 

t 

« 

MHABMF 
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INPTO 


Name;  INPTO 

Function:  input  for  old  job 


RESTRT 

DEBUGI 

NROTOR 

ANTYPE(3) 

OPREAD(lO) 

DEBUG (25) 

NPRNTI 


integer  parameter;  debug  print  control 


CASECM 

TMDATA 


I 


IWPTAl 


Name:  INPTAl 

Function:  read  airfoil  table  file 


DEBUG 

TITLE(20) 

IDENT(4) 

NMAX 

NAB 

NA(20) 

A(20) 

NMB 

NM(20) 

M(20) 

NRB 

R(ll) 

CLT(5000) 

GDT(5000} 

CMT(5000) 


THDATA 

AITABL 
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INPTRl 


Neime;  INPTRl 

F’unction:  read  rotor  anmelist 

DEBUG 

TITLE( 20) 

« 

t 

TWISTl(5l) 


TMDATA 

RIDATA 
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INPTWl 


Name  j  INPTWl 

Function!  read  wake  namelist 

DEBUG 

FACTWU 

« 

KWGS0(4) 

KFWG 

DQWG(2) 


TMDATA 

WIDATA 


GIDATA 
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INPTB 


Name :  INPTB 

Function:  read  body  namelist 

DEBUG  TODATA 

TITLE(20)  BDDATA 

DOFSYM(IO) 

lftaw  BADATA 

OPTINT 

ENGPOS  ENDATA 

« 

KEDAMP 
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INPTLl 


Name :  INPTLl 

Function:  read  loads  namelist 

DEBUG 

MHARML 

0PNCIS(4) 

MVIB 

ZETAV(3,i0,10) 


TMDATA 
LI DATA 

LA DATA 
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INPTF 


Name :  INPTF 

Function:  read  flutter  namelist  for  new  job 

DEBUG 

OPFLOW 

NAMEXR(3) 


IWDATA 

FLDATA 


ims 

Name :  INPTS 

Function:  read  flight  dynamics  namelist  for  new  job 

DEBUG  TMDATA 

NPRNTP  STDATA 

D0FPLT(2l) 


OPTRAN 


GCDATA 


INPTT 


Name  I  INPTT 

Function:  read  transient  namelist  for  new  job 

DEBUG  TMDATA 

NPRNTT  TNDATA 

OPLHDA 

OPTRAN  GCDATA 

« 

CHAG(5) 


-70- 


infix; 


Name :  INPTG 

Function:  read  flutter  namelist  for  old  job 
DEBUG 

ANTYPE(4) 

« 

NAHEXR(3) 


TODATA 

FLDATA 


-?1 


INPTU 


Name :  INITU 

Function:  read  flight  dynamics  namelist  for  old  job 

DEBUG  THDATA 

0PPRNT(4)  STDATA 

DCFPLT(2l) 

OPTRAN  GCDATA 

CMAG(5) 
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INPTV 


Name:  INFr\' 

'.  unction:  read  transient  namelist  for  old  jot 


DEBUG 

NP'^NTT 

NTONTP 

NPRN^^ 

NRS'URT 

TMAX 


TMDATA 

TNDATA 
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FILEI 


Name:  FILEl(mLE,RDWRT) 

Function:  read  or  write  input  file 

NFILE  file  unit  number 

REWRT  integer  parameter;  0  to  read  file, 

1  to  write  file 

TITLE D( 20) 

TrfLRl(20) 

TITLR2(20) 

BDDATA 

RiDATA 

TITLCS(20) 

R2DATA 

FILEID(4) 

TWDATA 

all 

all 

TMDATA 

all 

BDDATA 

all 

BA  DATA 

all 

SNDATA 

all 

LA  DATA 

all 

GCDATA 

all 

TNDATA 

all 

STDATA 

all 

FLDATA 

all 

RIDATA 

all 

WIDATA 

all 

GIDATA 

all 

LI  DATA 

all 

R2DATA 

all 

W2DATA 

all 

G2DATA 

L2DATA 

FILEJ 


Name: 

FILEJ( NFILE, RDWRT) 

Function:  read  or  write  trim  data  file 

NFILE 

file  unit  number 

RIWRT 

integer  paorameter:  0  to  read  file, 

1  to  write  file 

HPSI 

LEVEL! 

LEVEL?. 

TMDATA 

KNWl 

WIDATA 

MRGl 

MRLl 

KFWGl 

GIDATA 

KNW2 

W2DATA 

MRG2 

MRL2 

KFWG2 

G2DATA 

all 

TRIMCM 

all 

RODYCM 

all 

ENGNCM 

all 

GUSTCM 

all 

CONTCM 

all 

COm'CM 

all 

MNSCM 

all 

QBDCM 

all 

RTPICM 

all 

RHIGM 

all 

MDICM 

all 

INCICM 

all 

WKVICM 

all 

MNHICM 

all 

AESICM 

all 

MNRICH 

^  n  *1 

ATTI?** 

1V.4  1 

all 

QRICM 

all 

RTR2GM 

all 

RH2CM 

all 

MD2CM 

all 

INC2CM 

all 

WKV2CH 

all 

MNH2GM 

all 

ASS2GM 

all 

MNR2GM 

all 

AEF2GM 

all 

QR2GM 

-75- 


FILER 


Name:  FILER  (ri:WRT) 

Functions  read  or  write  restart  file 


Restart  file  structure; 


1)  case  header  record 

2)  input,  trim,  airfoil  data 

3)  task  header  record  —  ID,NREC 
(ID  =  2  for  flutter,  3  for  flight 
dynamics,  4  for  transient) 

4)  task  data  (NRBC  records) 

5)  repeat  ^3  and  ;^4  as  necessary 

6)  end  record  —  ID  =  0,  NREC  =  0 


RD^/RT 

integer  parameter:  0  to  read  file,  1  to  write  file 

RESTRT 

TITLGS(20) 

CASECM 

FILE1D(4) 

NROTOR 

con: 

TMDATA 

IDENT(4) 

TRIMCM 

TITLRi(20) 

TITLR2(20) 

TITLBD(20; 

RIDATA 

R2DATA 

TITLA1(20) 

BDDATA 

AF1ID(4) 

NMAXl 

cm  (5000 

) 

AITABL 

CDT1(5000’ 

GMTl(5000; 

) 

t 

TITLA2( 20) 

AF2ID(4) 

NMAX2 

A2TABL 

CLT2(5000) 

CDT2(5000) 

CMT2(5000) 

1 
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FILEF 


Name: 

FILEF(RDMRT) 

Function:  read  or  write  flutter  restart  file 

RMRT 

integer  parameter:  0  to  read  file,  1  to  write  file 

NRCTOR 

TMDATA 

OPFDAN 

FLDATA 

NBMl 

RTRICH 

NTWl 

NGMl 

NBM2 

RTR2CM 

Nma 

NGM2 

all 

FLMCM 

all 

STDGM 

all 

STMCM 

all 

MDICM 

all 

MD2CM 

all 

STRATA 

all 

GCDATA 
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FILES 


Name:  FILES (RDWRT) 

Function:  read  or  write  flight  dynamics  restart  file 


RDWRT  integer  parameter:  0  to  read  file,  1  to  write  fi]e 

STDGM 
STMGM 
STDATA 
GCDATA 


all 

all 

all 

cll 
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FILET 


Name:  FILET(RDWRT,ENDREC) 

Function:  read  or  write  transient  restart  file 


RDWRT  integer  parameter:  0  to  read  file,  1  to  write  file 

ENDREC  integer  parameter:  0  if  at  start  of  transient  record, 

1  if  at  end  of  record  (required  for  file  write  only) 

IT  WORK 

YN(7) 

DYN(7) 

DDYN(7) 

MTRAGE 
TRACE (14377) 

LEVELl  TMDATA 

LCVEL2 


all 

all 

all 

all 

all 

all 


TRANCM 
TNDATA 
GCIATA 
LI DATA 
L2DATA 
LADATA 
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FILEE 


Name;  FILEE(KEY) 

Function:  write  eigenvalue  file 


KEY 


integer  parameter  defining  case 


1 

2 

3 

4 

5 

6 

?-18 


start  file 
flutter,  const,  coeff.  (FLUTL) 
complete 
symmetric 
antisymmetric 

flutter,  periodic  coeff.  (FLUT) 
complete 
symmetric 
antisymmetric 
flight  dynamics  (STABL) 

6+IEQ  (IEQ  =  equation  type) 


TASK 

JGASE 

IDENT(4) 

CODE 

lamda(6o) 

MX2 

lmdap(6o) 

LMDAGP(60) 

MX2P 


CASECM 

TRIMCM 

\( constant  coefficients) 

TMDATA 

EIGVC 

\  (periodic  coefficients) 

\  (periodic  coefficients) 

EIGVP 
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TMDATA 


INITA 


Name:  INITA 

Function:  initialize  environment  parameters 
General  reference:  section  2.5 


OPUNIT 

ALTMSL 

TE3MP 

D2NSEI 

OPDENS 

DENSE 
A  LTD 
DRATIO 
CSOUND 


TMDATA 


TRIMCM 
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r  . 


MPSI 

MIIARM(2) 

MHARMF(2) 

OPTRIM 

OFGOVT 

LEVEL2 

D0F(54) 

IX)FT(8) 

VKTS 

VEL 

VTIP 

RPM 

COLL 

LATGYC 

LNGCYC 

PEDAL 

A PITCH 

AROLL 

ACLIMB 

A  YAW 

RTURN 

NROTOR 

XTRIM 

CXTRIM 

THETFT 

PIIIFT 

THETFP 

PSIFP 

THETAT 

PSIT 

IV  mrvrw/  ^ 
U  >  JJV  ux.  \  V/ 

DOMEGA 

DDZF 

VPILOT(5) 

TGOVRl 

TGOVR2 

NBLDl 

VTIPN 

RADIUS 

SIGMA 

GAMHAO 


CCNTCM 


RIDATA 


.  4 


•4 

sii 

li 
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o 


INITC 


OMEGA 1 

RTRICM 

OMEGA 2 

RTR2CM 

HMASS 

BODYCM 

TRATIO 

BDDATA 

CONFIG 

WEIGHT 

NBLD? 

R2DATA 

DRATIO 

TRTMCM 

DENSE 

GRAY 

C (TARG 

OFRTR2 

DPS  I 

FSCALE 

RSCALE 

NSCALE 

ISGALE 

GSCALE 

SSCALE 

CSGALE 

C0SPSI(36) 

SINPSI(36) 

KEPSI(21,36) 
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INITRl 


Name:  INITRl 

Function:  initialize  rotor  parameters 

Normalization  parameters:  section  2.6 

Aerodynamic  r,  ^^r;  section  2.4.1 

Tip  loss  factor:  section  2.4.5 

Linear  twist:  section  2.3.5 

Control  system  damping:  section  5.1.3 

Gimbal/teeter  spring  and  damping:  sections  2.2.12,  2.2  13 

Lag  damper:  section  2.2.16 

DEBUG 

MPSI 

rotor  degrees  of  freedom 

D0FT(4) 

LEVEL 

TRATIO 

DENSE 

CSOUND 

DRATIO 

QRTR(6) 

FHUBH(6) 

CLS 

CXS 

CTS 

CYS 

CPS 

CT 

GMX 

CMY 

BO 

BC 

BS 

K2T 

K2SI 

K2S0 

GAMMA 0 

SIGMA 

NBLADE 

RADIUS 

VTIPN 

TDAMPO 

TDAMPC 

TDAMPR 


TMDATA 


BDDATA 

TRIMCM 

QRICM 


WGICM 

RIDATA 
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INITRl 


NUGCO 

NUGSO 

GDAMPC 

G  DAMPS 

LDAMPC 

LDAMPM 

L  DAM  PR 

MRS 

MRM 

RAE(3l) 

MRA 

BTIP 

CPTIP 

TWISTA(30 

TWISTI(51 

RI(51) 

MR  I 

INFL0W(6) 

LINTW 

TWISTL 

OMEGA 

GLAG 

HLD 

DZLD 

GGS 

CGC 

NUGC 

NUGS 

CTO 

CTC 

CTR 

MTIP 

GAMMA 

CMEAN 


RIDATA 


RTRICM 


NBM 

NTM 

NGM 

NBMT 

RA(30) 

DRA(30) 

F11P(30) 

CTOLD 

GMXCLD 

CMYOLD 

VIND(3.30,36) 

LAMBDA 


WKVICM 
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INIRl 


VINT(3,30,36)  WKVICM 

VORH(3,36) 

LAHBDI 

vwb(3,36) 

vht(3,36) 

vvt(3,36) 

voff(3,36) 

LAHBEW(3) 

LAMBDH(3) 

LAIffiDv(3) 

LAMBIX)(3) 

EINTW(3) 

EINTH(3) 

EINTV(3) 

STATE(30,36,3)  AESICH 

dclmaxOo,36) 

dcemax(30,36) 

DCMMAX(30,36) 

alpha(3o,36) 

BETA(21,10)  MNRICM 

THETA(21,5) 

BETAG(2l) 

PHI(10,16) 

PSID(10,6) 

QSSTAT(lO)  MNSCM 

PISTAT 

PESTAT 

FORCE( 16,36)  AEFICM 

FHUB(6,36j 

TORQUE(36) 

T750LD  MDICM 

NBMOLD 

NTMOLD 

VGUST(3,30,36)  GUSTCM 

'  GUSTH(3) 
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o 


o 


INITB 


Name:  INITB 


Function;  initialize  airframe  parameters 


Position  of  aircraft  components; 
Rotation  matrix  Rgp:  section  4.1.2 


r,  Rgp  without 


section  4.1.5 
sections  4.1.3, 


4.1.5 


(for  wind  tunnel  trim  case) 

Control  matrix  Tpp„;  section  4.1.6 
Aircraft  inertia:  section  4.2.4 
Airframe  elastic  modes; 

a)  pitch/mast-bending  coupling  (KMST):  section  4.2.3 

b)  mode  shape  at  hub  (AHODE):  section  4.2.2 

c)  mass,  sprijig,  damping:  section  4.2.4 

d)  aerodynamic  damping  and  control:  section  4.2.? 


Initialization  (for  wind  runnel  case) 

"fv  - “  '‘fE  ■ 


«e-"f 

-  H  (  V  X  )  Rg  =  -  M  V  (tp  X  ) 
("if^)  Rg  kp  =  -  V  jj; 


G  =  -  M*  g  (kp  X  ) 


DEBUG 

VEL 

D0F(16) 

airframe  degrees  of  freedom 

TMDATA 

GRAY 

GAMMA 

SIGMA 

IB 

OMEGA 

NBLADE 

RADIUS 

reference  rotor 

TRIMCH 

P21MR1 

0. 

RTRICM 

P21MR1 

(rad) 

P21MR2 

(rad) 

RTR2CM 

P21WR2 

(rad) 

ROTATl 

0PHVB1(3) 

RIDATA 
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C 


INITB 


R0TAT2 

0PHVB2(3) 

R2nATA 

VGVrBV(3)  gust  in  velocity  axes 

VGHTV(3) 

GUSTCH 

VGVTV( 

:3) 

QWB(6; 

qht(6; 

QVt(6, 

) 

) 

) 

QBDCM 

AM0DE1(6,10) 

• 

• 

BODYCM 

• 

VSIDE 

TITLEUO) 

* 

• 

• 

)i)FSYM(lO) 

BDDATA 

DRGIW 

BA  DATA 
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INITE 


Name;  INITE 

Function:  initialize  drive  train  parameters 

Engine  inertia  and  control;  sections  4.3.1,  4.3.2 
Governor  parameters  (dimensionless);  section  4.3.3 
Drive  train  spring  constants:  section  4.3.2 


DEBUG 

OPENGN 

D0F(6) 

drive  train  degrees  of  freedom 

TMDATA 

TRATIO 

BDDATA 

NBLADE 

IB 

OMEGA 

reference  rotor 

TRIMCM 

ENGPCS 

EVDATA 

THRTLC 

I  ENG 

KMASTl 

KMAST2 

KICS 

KENG 

KPE 

KPl 

KP2 

TIE 

Til 

T12 

T2E 

T21 

T22 

QTHRTL  ENGNCM 

lENGS 

KMTl 

KMI2 

KMR 

KMEl 

KME2 

KPGOVE 

KPGOVl 

KPGOi^2 

TIGOVE 

TlGOVl 

T1GOV2 

T2GOVE 

T2G0V1 

T2GOV2 

NDM 
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Function;  check  for  fatal  errors 

MPSI 

LSV2L 

NBLADE 

MPA 

RAE(31) 

MHI 

RI(51) 

RRCOT 
INFLOW (6) 

MRG 

NG(30) 

HRL 

NL(30) 

KNW 

RA(30) 


TMDATA 

RIDATA 

WIDATA 


MRTO 


other  rotor 


RTRICM 

R2DATA 


PRNTJ 


Name :  PRNTJ 

Function:  print  job  input  data 

TMDATA 
CASECM 
UNITNC 


FILEID(4) 

all 

all 
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PRNTC 


Name :  PRNTC 

Function:  print  case  input  da,ta 


JCASE 

JOB 

START 

FILEI  i(4) 
TITLCS( 20) 
CODE 

ANTYPE('3) 

CPUNIT 

cm  IN 

NRC  TOR 

VKTS 

VEL 

RPH 

''TIP 

AmiSL 

TSI'iP 

CPGRND 

HAGL 

A  FLAP 

OPENGN 

OPGOVT 

RTiJRN 

LE'/BLl 

LS'''SL2 

1X'F(54) 

DOFT(R) 

MI'S  I 
MHARM 
MHi\RMF 
OPDENS 

IDENT(4) 

DENSE 

DRATIO 

CSOUND 

ALTD 

TITLBD( 20) 

WEIGHT 

FSCG 

WLGG 

BLCG 

CONFIG 

ATILT 


CASECH 


TI'IDATA 


TRIMCM 


BDDATA 
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qr<-£wiW*' 


PRNIG 


CWS 

NDM 

TITLAl(20) 

AF1ID(4) 

TITLA2(20) 

AF2ID(4) 

'  TITLRl(20) 
TYPEl 
RADUSl 
NBLDl 
SIGMAl 
INFLWi(6) 

0PHVB1(3) 

OPSTLl 

OPYAWl 

OPCMPl 

OPUSLl 

ROTATl 

HINGEl 

ELAGl 

EFLAPl 

GAMMA 1 

OMEGAl 

MTIPl 

CMEANl 

IBl 

NBMl 

NTMl 

NGMl 

NBMTl 

TITLR2(20) 


BODYCM 

ENGNCH 

AITABL 

A2TABL 

RIDATA 


HTRICM 


RPDATA 


EFLAP2 
GAMMA 2 


RTR2CM 


NBHT2 
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Function:  print  trim  input  data 
FILEID(4) 

HHARMF( 2) 


TODATA 


PRNTOl 


Name ;  PRNTRl 

Function:  print  rotor  input  data 


NBM  RTRICM 

NTH 

NGH 

RA(30) 

DRA(30) 

1-TIPbo) 

TITLE(20)  RIDATA 

TCISTI(51) 
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PRNTWl 


Name :  PRNTVI 1 

Function!  print  wake  input  data 

MPSI 

TCDATA 

LEVEL 

FACTC'R 

« 

WIDATA 

• 

K’rfGS0(4) 

KFVJG 

GIDATA 

» 

DQWG(2) 
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PRNTB 


Nar>e:  PRNTB 

Function:  print  body  input  data 

NROTC'R 

TITLE( 20) 

* 

IX)FSYK(lO) 

LFTAW 


TMDATA 

3DDATA 


BA DATA 


OPTINT 

ENGFUS 


ENDATA 


KEDAHP 


PRNTF 


Name:  PRNTF 

Function:  print  flutter  input  data 

IDENT(4) 

CONFIG 

CPFLCW 

< 

GPUSLD 


TRIHCH 

BDmTA 

FLDATA 


PR  NTS 


Name :  PRNTS 

Function:  print  flight  dynamics  input  data 

IDENT^4) 

NPRNTP 

CHS(3) 


TRIMCM 

STDATA 


PRNTi 


t  -X.--  -A-''  •TV^V’. 


Name :  PRNrr 

Function:  print  transient  input  data 

IDENT(4)  TRIMCM 

NPliNTT  TNDATA 

(^IWDA 
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PRNTG 


Name:  PRNTG 

Function:  print  transient  gust  and  control  input  data 

NROTDR 

OPTRAN 


TMDATA 

GGDATA 
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TRIM 


Neune:  TRIM 
function:  trin 

General  reference:  sections  5*3«5f  5*3*1 


RESTRT 

RSWRT 

CPRTR2 

LS'.'SLl 

LSVEL2 

ITERU 

ITERR 

TTERF 

NFRNTT 

NPRNTP 

NPRNTL 


GASEGM 

raiMcn 

Tl'i  Al'A 
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TRIMI 


Name:  TRIMI ( LEVELl ,LEVEL2) 

Function:  calculate  trim  solution  by  iteration 
General  reference:  section  5*3 •! 

Codes : 

control  number  (C)  =  123  ^••5  6  7  8 

control  =  Sq  4  Sp  ©PT  ‘fpT  %p  ^pp 

test  number  (t)=1  23^-  36?  8  9  10 

test  =  none  F  M  F^F^^  Gp  C^,  0^0, 


OPTRIM 

C(i) 

T(i) 

1  2  3  5  6 

211311 

1  2  3  5  7 

211311 

3 

7 

123^568 

2113116 

4 

7 

123^578 

2  113  116 

5 

1  3  5 

4  1  5 

6 

MM 

1  3  5  8 

4156 

7 

R 

HI 

9 

HI 

10 

HI 

11 

mM 

1 

7 

12 

1 

9 

7 

13 

1 

1 

6 

14 

2 

2  3 

8  9 

15 

3 

1  2  3 

7  8  9 

16 

3 

123 

11  1  1 

17 

3 

1  2  9 

11  1  1 

18 

1239 

10  1  8  9 

19 

1  2  3 

11  1  1 

20 

3 

1  2  9 

11  1  1 

21 

4 

1239 

10  1  8  9 

22 

1 

3 

8 

23 

1  3 

7  8 

24 

1  3 

10  1 

25 

2 

1  9 

10  1 

26 

3 

1  3  9 

10  1  ® 

27 

2 

1  3 

10  1 

28 

2 

1  9 

10  1 

29 

3 

1  3  9 

10  1  8 

0' 


LEVELl 

LEVEL2 


waike  analysis  for  rotor  "1  and  rotor  "2: 
0  for  uniform  inflow,  1  for  prescribed 
wake,  2  for  free  wake 


TRIMI 


DEBUG 

TMDATA 

CPTRIM 

CTTRIM 

CYTRIM 

BSTRIM 

BCTRIM 

OPTRIM 

MTRIM 

HTRIMD 

FACTOR 

ITERM 

ITERC 

DELTA 

EPTRIM 

OPGOVT 

CXTARG 

TRIMCH 

GRAY 

COUNTT 

CNTRLZ(ll) 

BDDATA 

CWS 

BODYCM 

KE(3) 

VXREKF(3) 

TCFE(ll,5) 

COUNTl-I 

CO  NY  CM 

GOUNTG 

NBLDl 

RIDATA 

ROTATE 

NBLD2 

R2DATA 

GAMMA 1 

RTRICM 

OMEGA 1 

IBl 

GAMMA 2 

RTR2CM 

OMEGA 2 

IB2 

VCNTRL(il) 

CONTCM 

THETFT 

PHI  FT 

THETFP 

PSIFP 

THETAT 

DPSIF  '^p 

VPILOT(5) 

TGOVRl 

TGOVR2 
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TRIMI 


QRTR1(6) 

GLS 

GXS 

GTS 

GYS 

GFS 

BETAG 

BETAS 

GQSl 

V-  = 

G^tr- 

QBIGM 

QRTR2(6) 

GQSR 

\J<^ 

QRGGf’ 

QWB(6) 

qht(6) 

QVT(6) 

QBDGM 

raiMP 


Name :  TRIMP( LEVELl , LEVEL2, ITER , ITERH) 
Function:  print  trim  solution 


LEVELl 

LEVELS 


ITER 

ITERM 

GPTRIH 

CTTRIM 

CYTRTM 

BOTRIM 

BSTRIM 

OPTRTM 

HTRIM 

E/TRIM 

CPGC'/T 

CCLL 

LATCYC 

LNGCYG 

PEDAL 

A PITCH 

A  YAW 

A  ROLL 

A CLIMB 

CXTARG 

GRAV 

COUNTT 

CPRTR2 

NBLDl 

TYPEl 

NBLD2 

TYPE2 

GAMMAl 

OMEGAl 

IBl 

GAMMA 2 
OMEGA 2 
IB  2 

CWS 

KE(3) 

VXREKP(3) 


wake  analysis  for  rotor  and  rotor 
0  for  uniform  inflow,  1  for  prescribed 
wake,  2  for  free  wake 

iteration  number 

maximum  number  of  iterations 


THDATA 


TRIMCH 

RIDATA 

RSDATA 

RTRICM 

RTR2CM 

BODYCM 
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TRIMP 


VCNTRL(ll) 

THE TFT 

PHIFT 

THETFP 

PSIFP 

THETAT 

PS  IT 

• 

GCNTCM 

D)’SIF 

VPILC'T(5) 

TGCVRl 

TGOVR2 

QRTR1(6) 

GLS 

CXS 

GTS 

GYS 

GPS 

BETAG 

BETAS 

GqSl 

=  Cp/<7- 

QRICM 

QRTR2(6) 

CQS2 

II 

o 

QR2CM 

CqA- 

QWB(6) 

QBDCM 

QHT(6) 

QVT( 6) 

-lOR- 


FLUT 


Name :  FLUT 
Function:  flutter 

General  reference:  sections  5 *3 *51  5 *3 *6 


.'^SWRT 

RESTRT 

0PRTR2 

NBLADE 

CPFLOW 

GKYMM 

CPFDAN 

MPSIPG 

NINTPC 

NBLDFL 

A2(6'+00) 

MGA 

MXFD 


CASECM 

TRIMCM 

FLDATA 

FLMCM 

STMCM 
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FLU-TH 


Names  FLUn^PSl) 

Functions  calculate  flutter  matrices 
General  references  section  6.3.1 


Inflow  dynamics s 


sections  6.1.5,  2.4.3 

DLDT=  z7 


DLDM  = 


Z-V 


TT  =  Tl-r 
TM  =  *x:y^ 

DLDZ  =  Si 

ZK  =  tv 

Drive  train  equations s  section  6.2.3 
Construct  flight  dynamics  matrices s  section  5 *3 *3  also 
(only  if  rigid  body  degrees  of  freedom  present) 
Symmetric/antisymmetric  matrices s  section  6.3.3 


PSI  ‘V*  (for  periodic  coefficients) 

DEBUG 

OPENGN 

0PRTR2 

DOFSYM(IO) 

TRATIO 

CONFIG 

NEM 

REULER(3,3) 

KE(3) 

RKUB1(3) 

RHUB2(3) 

AM0DE1(6,10) 

AMODE2(6,1o) 

KMSTC1(10} 

KMSTSl(lO) 

KMSTC2( 10} 

KMSTS2(lO) 

MVXRE(3,3) 

TCFE(11,5) 

KIGOVE 

KIGOVl 

KIGOV2 


THDATA 

TRIMCH 

BDDATA 


BODYCM 


ENDATA 
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FLUTO 


GSB 

GSl 

QTHHTL 

I  ENG 

QBDAMP 

KMIl 

KNI2 

Krn 

KMEl 

KMS2 

KPGC’^E 

KPGe'’l 

KPGCV2 

TIGC' E 

TIGC’  1 

VlGCv-2 

T2GC’VE 

T2Grvi 

T2GCV2 

MENG22 

HSNG33 

SENG22 

SENG33 

PA DUS 1 

NBLDl 

KFLNDl 

KHLMDl 

SIGKAl 

FXLHDl 

fylmdi 

KINTHl 

KINTPl 

FMLMDl 

OMEGA 1 

NTMl 

NBMl 

NGMl 

MUXl 

MUYl 

HUZl 

GAMMA 1 

IBl 

RGUST1(3,3) 

CHUB1(6,16) 

CBHUB1(3,3) 

CHUBT1(16,6) 


BNDATA 

ENGNCM 


RIDATA 


RTRICM 
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RADUS2 


FMLMD2 
OMEGA 2 


GHUBT2(i6,6) 

KPBl(lO) 

KPGl 

KPB2(10) 

KPG2 

TlCl 

TlSl 

T1C2 

T1S2 

LAMBDI 

COSEl 

ZAGLl 

LAMBD2 

COSE2 

LAMBD2 

CTSl 

CTS2 

DERIV(7,21) 

DRVR1(7,21) 

DRVWB(7,21} 

DRVHT(7,21) 

DRVVT(7,21) 

A2FD(7.7) 


20^7^  a 

y  20^^- a 


MXFD 

CPFLOW 

OPSYHM 

NBLADE 

OPSAS 

KCSAS 

KSSAS 

TCSAS 

TSSAS 

0PTCRS(2) 

CPGRND 

KASGB 


FLUTM 

R2DATA 


RTR2CM 


MDICM 


HD2CH 


GONTCH 


WKVICH 

WKV2CM 

QRICM 

QR2CH 

STDCH 


STMCM 


FLDATA 
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FLUTM 


I.  ij.'' 


D0F(80) 

CON( 26) 
GUS(3) 

A 2(6400) 

HGa 

A2A(16,16) 

« 

« 

BLA(16,2) 

A2R1(30,30) 

DGR1(6,3) 

A2R2(30,30) 

DGR2(6.3) 


FLDATA 

FLMCH 

FLMACH 

FLMICM 

FLM2CM 
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FLUTB 


Name :  FLUTB 

Function:  calculate  flutter  aircraft  matrices 
General  reference:  section  6.2.2 


0PRTR2 

TRIMCM 

NSH 

BDDATA 

IBODY(3.3) 

BCDYCM 

HSTAR 

MVXRE(3,3) 

GHTRX(3,3) 

rfv(3.3) 

AMASS(IO) 

AIAMPS(IO) 

ASPRNG(IO) 

ADAMPA(lO) 

ACNTRL(4,10) 

DELTA 

('PRINT 

rLlMTA 

D'.'BODY(6) 

DDZF 

CNTRL(4) 

(Sf 

''A  NT'F 

GWB(3) 

gust  in  F  axes 

GUSTCF 

GHT(3) 

GV'?(3) 

QWB(6) 

QHT(6) 

OBDC."! 

QVT(6) 

A2(16,16) 

FLMAT’ 

DRVWB(7,2l) 

UR' KT(7,21) 

DR''Vt(7,21) 

STDGM 

LMDAW1(3) 

LMDAHl(3i 

LMDA”1(3 

EINTWi(3J 

EINTHi(3] 

EINT\'1(3> 

WKVICM 

LMDAW2(3) 

WKV2CM 

EINTV2(3) 
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KLUTRl 


Name:  ii'LUTRl(PSl) 

Function:  calculate  flutter  rotor  matrices 
General  reference:  sections  6.1.6,  6.4 

Azimuthal  s"jiimations : 


to 

tsi 

VW  K  J 

at 

.  1,  Ztr- 

=  7;r 

for  periodic  coefficients 

J- 

"T 

ir 

at 

4^. 

=  i  — 

<3  TT 

for  constant  coefficient 

) 

approximation 

(section  6.1.?) 

Reorder  hub  reactions:  _A.^  equation  multiplied  by  2  to  get  (  -’S'ZC^v'a) 
Inflow  dynamics  due  to  velocity  perturbations:  sections  6.1.4,  6.1-6 

PSI  (periodic  coefficients  only) 

CPFLCW  FLDATA 

NPSICG 

NBLDFL 

KBM  RTRICM 

KTO 

NGH 

GAMMA 

Nl'GC 

NIXIS 

CGC 

CGS 

CTO 

CTC 

GTC 

MUX 

MUY 

MUZ 

NBLD  RIDATA 

GSB(IO) 

GST(5) 

KMLMDA 

KFLMDA 

NU(IO)  •  MDIGM 

WT(5) 

WTO 

WTG 

WTR 

KPB(IO) 

KPG 
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FLUTRl 


LA^^BD^ 

WKVICM 

CTS  'J( 

QRIGM 

TIG 

CONTCM 

TIS 

A2(30,30) 

• 

FLMICM 

DG(6,3) 

MASSE 

• 

FLINCH 

• 

spq(5,io) 

MQU(IO) 

• 

FLAECM 

MPDP(5.5) 
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FLUTIl 


Name;  FLUTIl ( PST ) 

Function:  calculate  flutter  inertia  coefficients 
General  reference:  section  6.1.3 


PST  M-' 

DEBUG  TMDATA 

ECFT(4) 

GLAG  RTRICM 

K3M 

KT5-I 

NBMT 

BETA(21,10)  MNRICM 

ETAPH(2,10)  MDICH 

MB  INCICM 

• 

W(5,10,4) 

HASSBL  FLINCH 


SPQL(5,10) 


FLUTAl 


Namet  FLUTAl (PSi) 

Function:  calculate  flutter  aerodynamic  coefficients 
General  reference:  section  6,1.4 


Perturbation  section  forces:  without  c/c  factor 
Aerodynamic  coefficients:  FZO  =  ,  FXO  = 

PSI  4^ 

DEBUG 

D0FT(4) 

MPSI 

DPSI 

MRA 

GH0RD(30) 

XA(30) 

XAC(30) 

CPGOMP 

OPYAW 

CPSTLL 

RFA 

RA(30) 

DRA(30) 

CMEAN 

FTIP(30) 

N3MT 

KBH 

KTH 

MTIP 

MUX 

MUY 

MUZ 

ETA(2. 10,30)  bending  modes  at  r, ,  i  =  1  to  MRA 
ETAP(2,10,30) 

ETAPP(2,10,30) 

ZETA(5»30)  torsion  modes  at  r, ,  i  =  1  to  MRA 

ZETAP(5,30)  ^ 

DELI 

DEL2 

DEL3 

DEL4 

DEL5 

DALPHA 

DMACH 

OPUSLD 


TMDATA 


TRIHCM 

RIDATA 


RTRICM 


MDICM 


FLDATA 
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FLUTAl 


BETA(21,10)  MNRICM 

DCl^(30,36}  AESICM 

DCDI]S(30,36) 

DCMI]S(30,36) 

SAVE(30,36,19) 

XAPQ(2,5.4,30)  INCICM 

MQU(IO)  FLAECM 

• 


NPDP(5,5) 
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KLUTL 


Name :  FLUTL(  ID,  A2,A1  ,A0  ,B  .'IX  ,MX1  ,HV  ,MG  ,  DOFl ,  NAMEX  ,NAKEV) 
Function:  analyze  flutter  constant  coefficient  linear  equations 

Vibration  point  location:  sections  4.1.3,  4.1.5 


ID 


A2(MX*M.<) 

Al(MX*MXj 

A0(MX*HX) 

B(MX*MV) 

MX 

MXl 

MV 

MG 

DOFl (MX) 

NA’!EX(MX) 

NAHEV(HV) 


■EL‘-(3) 

GHA'’ 
r  MEGA 
RADIUS 

!<SGG 

BLCG 

WLCG 

NSM 


THE'H-T 
PHI  FT 
THETA T 
PS  IT 

ANTYPE(4) 

NAMEXR(3) 


problem  identification:  1  for  complete  dynamics, 
2  for  symmetric,  3  for  antisymmetric 

coefficient  matrices 


control  matrix 

number  of  degrees  of  freedom 

number  of  first  order  degrees  of  freedom 

number  of  controls 

number  of  gust  components 

integer  vector  designating  first  order  degrees 
of  freedom 

vector  of  variable  names 
vector  of  control  names 

reference  rotor 
reference  rotor 


BODYG" 

TRIMCH 

BDDATA 

CONTCM 

FLDATA 
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STAB 


Name :  STAB 

Function;  flight  dynamics 

General  reference:  sections  5*3*5i  5*3*3 

RESTRT 

RSWRT 


CASECM 
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^v.  -  t.»^--i.-<igatRajiagaRTO;ggMqgwyi 


STABH 


Name:  STABH 

Function:  calculate  flight  dynajnics  stability  derivatives  and  matricec 
General  reference:  section  5 *3 *3 

Print  during  stability  derivative  calculations: 

a)  increment:  1st  number  dimensionless,  2nd  number  dimensional 

b)  motion  and  controls:  1st  number  dimensionless,  2nd  number 

dimensional 

1^  angular  velocity  =  deg/sec 

2)  linear  velocity,  gust  velocity  =  ft/sec  or  m/sec 

3) 

2  2 

4)  2p  =  ft/ sec  or  m/sec‘^ 

5)  controls  =  deg 

c)  generalized  forces:  moments  and  forces  in  ?i2C/c— a  for:.i 

(rotor  #1  paxameters,  body  axes); 
torque  in  0^/4- a  form  (rotor  n^l 

parameters)  ^ 


MPSI 

TMDATA 

LEVELl 

LEVEL2 

DEBUG 

0FRTR2 

THIHCM 

LSCALE 

FSCALE 

NBLDl 

RIDATA 

MRAl 

TYPEl 

IBl 

RTOICM 

CHUB1(6,16) 

CHUBT1(16,6) 

OMEGA 1 

NBLD2 

R2DATA 

MRA2 

TYPE  2 

IB  2 

RTR2CM 

CHUB2(6,16) 

GHUBT2(16,6) 

OMEGA 2 

ibody(3,3) 

bodycm 

MSTAR 

MVXRE(3,3} 

GMTRX(3,3) 

TCFE(ll,5) 
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STABM 


CONFIG 

BDDATA 

qRTRl(6) 

QRICM 

CQSl  -6  20^^  a 

QRTR2(6) 

QR2CM 

CQS2 

101 

INClCtl 

102 

INC2CM 

IRSTAR 

ENGNCM 

QTHRTL 

QEDAMP 

KPGOVE 

KFGOVl 

KPG0''2 

KIGOVE 

ENDATA 

K ICO VI 

KIGCV2 

NPRNTP 

STDATA 

N^RNTL 

I'lERS 

OPLMEA 

DELTA 

D0F(7) 

CCN(l6) 

GUS(3) 

VGWBV(3) 

GUSTCH 

VGHTV(3) 

VGVTV(3) 

vgrtriO,30,36) 

vgrtH2(3,30,36) 

VGHUBlOj 

VGHL'B2(3) 

/CNTRUllJ 

11  i  on 

r/ 'BODY  (6) 

r>'.-'3GA 

I^jZ" 

qwb(6} 

QBDCM 

QHT(65 

QVT(6) 

DERIV(7,2l) 

STDCH 

• 

• 

DRVVT(7,21) 

A2FD(7,7) 

STMCM 

M7JD 
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STABD 


Name:  STABD 

Function:  print  stability  derivatives 
General  reference:  section  5*3*3 

Options:  a)  rotor  coefficient  form,  =  ^f'2C/«r‘a 

stability  derivative  form,  X  (acceleration) 
c)  dimensionless  or  dimensional 

Dimensions : 

a)  force  or  moment 

forces (FF) 

M*X  form 
X  form 

b)  subscripts 

acceleration  { *z )  =  (fa) 

angular  velocity  =  n 

linear  velocity  =  JTIR  (FV) 

controls  =  57*3 

gust  velocity  =  riR  (FV) 

TASK 

D0FFD(7) 

GONFDf  16) 

GUSFB(3) 

NAMEV(19) 

ISTAR(3,3) 
nSTAR 

IRSTAR 

NBLADE  reference  rotor 

IB 

OMEGA 
RADIUS 

CPPRNT(4) 

DRVRl(y,2l) 

DRVR2(7,21) 

DRVWB(7,2l) 

DRVHT(7,21) 

DRVVT(7,21) 


CASECH 

STMCM 


BCDYCM 

ENGNCM 

TRIMCM 


STRATA 

STDGM 


STABE 


Name:  STABE 

Function:  calculate  flight  dynamics  equations 


DEBUG  TODATA 

CIIEGA  reference  rotor  TRIMCM 

EQTYPE(12)  STDATA 

NCSAS 

KSSAS 

TCSAS 

TSSAS 

A2FD{1^9)  SniC'-J 

MXFD 

CPSYMH  FLDATA 

C  KASF 

TASK  CASECh 
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STABL 


Name:  STABL(IEQ,A2,A1  .AO.B.MX.MXl.HV.HG.lOFi  .NAI-EX.NAHEV.DOF.CON) 
P'unction:  analyze  flight  dynamics  linear  equations 

Vibration  point  location:  sections  4.1.3,  4,1.5 
Numerical  integration  of  transient:  sections  5 *3 *2,  5 *3 *3 
(see  also  program  TRAN) 


lEQ 

equation  type  identifier 

A2(HX*MX) 

ai(mx^^hx) 

A0(HX*MX) 

coefficient  matrices 

b(hx^<mv) 

control  matrix 

MX 

ntunber  of  degrees  of  freedom 

MXl 

nujnber  of  first  order  degrees  of  freedom 

MV 

number  of  controls 

HG 

number  of  gust  components 

DOFl(MX) 

integer  vector  designating  first  order  degrees 
freedom 

of 

NAMEX(MX) 

vector  of  variable  names 

namev(hv) 

vector  of  control  names 

D0F(7) 

integer  vector  designating  degrees  of  freedom 

used 

C0N(19) 

integer  vector  designating  controls  used 

OMEGA 

RADIUS 

GRAV 

reference  rotor 

TRIMCM 

VELF(3) 

BCDYCM 

VGHUB1(3) 

VFrRAN(5) 

GUSTCM 

FSCG 

BDDATA 

WLCG 

BLCG 

THETFT 

PHIFT 

THETAT 

PSIT 

D'^0DY(6) 

DOMEGA 

CCNTCM 

NPRNTT 

STDATA 

D0FPLT( 21 ) 
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STABP 


Name :  STABP( TIM , IT . YN , DYN , DDYN , DOF ) 

Function:  print  flight  dynamics  transient  solution 
General  reference:  section  5 *3 *3 


Print  during  numerical  integration  (in  STABL) : 

a)  controls  in  deg 

b)  gust  velocity:  1st  number  dimensionless,  2nd  num'.Br 

dimensional 

c)  aircraft  motion:  1st  number  dimensionless,  2nd  number 

dimensional 


1)  displacement  =  deg,  ft  or  m 

2)  velocity  =  deg/sec,  ft/sec  or  m/sec 

3)  acceleration  =  deg/sec'^,  g 

4)  inertial  axes  =  deg/sec,  g 


AANG 


ALIN 


(V  X  ) 


TIM 

IT 

YN(?) 

DYN(7) 

umN(?) 

DCF(7) 

GRAY 

LSCALE 

FSCALE 


time  (dimensionless) 
time  count 


T  ^F 

(4>V 

( 'ip 


^P 


F  F 


V  “f 


p  ^  Yp  "F 


^  nr 

F  •'F  "F  'I 


integer  vector:  0  if  degree  of  freedom  not  used 


TRIMCM 


TSTEP  STDATA 

TMAX 

NPRNTT 
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STABP 


VGHUB1(3)  GUSTCM 

VPTRAN(5) 

MSTAR  BODYCH 

HVXRE(3,3) 

reulerO.3) 
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TRAN 


Ncime:  TRAN 

Function:  transient 

General  reference:  sections  5*3*5«  5 ‘3 *2 

RESTRT 

RSWRT 

CASECM 

LE'^ELl 

LE''EL2 

TI-IRATA 

DVBCDY(6) 

D'-'MEGA 

GONTCN 

MVXRE(3,3) 

MSTAR 

IBODY(3,3) 

BODY  CM 

OMEGA 

reference  rotor 

TRIMCM 

QRTR1(6) 

QRICH 

CQSl 

QRTR2(6) 

Hi 

-Y2CQ/c-a 

QR2CM 

CQS2 

QWB(6) 

QHT(6) 

QVT(6) 

QBDGM 

QTRIM(6) 

CQSTl 

CQST2 

IB0DYI(7,7) 

TRANCM 

NPRNTT 

NPRNTP 

NPRNTL 

NKSTRT 

TMAX 

TSTEP 

OPPLCT 

D0FPLT( 21 ) 
DOF(7) 

TNDATA 

101 

INGICM 

102 

INC2CM 

CHUB1(6,16) 
GHUBT1(16,6) 
OMEGA 1 

IBl 

RTRICM 
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TRAN 


CHUB2(6,16) 
CHUBT2(16,6) 
OMEGA 2 
IB  2 

NBLDl 

NBLD2 

IRSTAK 


RTR2CM 


RIDATA 

R2DATA 

ENGNGM 


-130- 


IRANI 


Name:  TRANl(y,DY,DDY) 

Function:  calculate  transient  acceleration  for  numerical  integration 
General  reference:  section  5 *3 *2 

Y(7)  (  'I’p  '^p  Xp  Yp  Zp 

DY(7)  (4>„  e„  %  X  y  z  ) 

DDY(7)  (  ‘4^  ep  iPp  xp  yp  Zp 


LBVELl 

LEVEL2 

DEBUG 


0PRTR2 

MVXRE(3,3) 

GMTRX(3,3j 

TGFE(11,5) 

GNmZ(ll) 

QTHRTL 

QEDAMP 

KPGOVE 

KKX)V1 

KPG0V2 

KIGOVE 

KIGOVl 

KIGOV2 

IBl 

OMEGAl 

NBLDl 

IB2 

0MEGA2 

NBLD2 

QRTR1(6) 

GQSl 

QRTR2(6) 

GQ32 

QWB(6) 

QHT(6j 

QVT(6) 


-  ^2G|^/v-a 

-  5  20^/^a. 


TMDATA 

TRIMGM 

BODYGM 

BDDATA 

ENGNGM 


ENDATA 

RTRIGM 

RIDATA 

RTR2GM 

R2DATA 

QRIGM 

QR2GM 

QBDGM 
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TRANI 


IX)F(7) 

OPSAS 

KCSAS 

KSSAS 

TCSAS 

TSSAS 

ITERT 

OPLMDA 

QTRIM(6) 

CQSTl 

CQST2 

IB0DYI(7,7) 

DCSAS 

DSSAS 

TTGOV 

TIGOV 

T2GOV 

V:'T^AN(5) 

VCNTRLfll) 

DVB0DY(6) 

DOMEGA 

DDZF 

VPILOT(5) 

TGOVRl 

TGOVR2 


TNDATA 


TRANCM 


GUSTCM 

CONTCM 
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TRANP 


Names  TRANP(TIM,IT,YN,DYN,DDYN) 
Function:  print  transient  solution 
General  references  section  5.3.2 


Print  notes: 

a)  controls  in  deg 

b)  gust  velocity  dimensionsal 

c)  aircraft  motions  1st  number  dimensionless,  2nd  number 

dimensional 

1)  displacement  =  deg,  ft  or  m 

2)  velocity  =  deg/sec,  ft/sec  or  m/sec 

3)  acceleration  =  deg/sec^,  g 

4)  inertial  axes  =  deg/sec,  g 

d)  generalized  forces:  moments  and  forces  in  »2C/^a  form 

(rotor  #1  parameters,  body  axes); 
torque  in  -  '?ic,  4r-a  form  (rotor  #i 
parameters)  ^ 


AANG  =  j 


ALIN  = 


^ody 


TIM 

time 

(dimensionless) 

IT 

time 

count 

YN(7) 

if. 

f  F  y  V  yp 

DYN(7) 

"^F  ^F  h  h 

DDYN(7) 

LBVELl 

LEVEL2 

FSCALE 

LSCALE 

GRAY 

0PRa'R2 

dp  %  ^  yp 

TMQATA 


TRIMCH 


TRANP 


ITERT 

OPLMDA 

TSTEP 

THAX 

MSTAR 

REULER(3.3) 

HVXRB(3,3) 

GMTRX(3,3) 

QTHRTL 

QEDAMP 

VGWBV(3} 

VGHTV(3) 

VGVTV(3) 

VGHUBlh) 

VGHUB2(35 

VPTRAN(5; 

NBLDl 

TYPEl 

331 

OMEGAl 

NBLD2 

TYPE2 

IB2 

OMEGA2 

QRTR1(6) 

CQSl 

QRTR2(6) 

CQS2 

QWB(6} 

QHT(6) 

QVT(6) 

VCNTRL(ll) 

VPILOT(5) 

TGOVRl 

TGOVR2 

QTRIM(6) 

CQSTl 

CQST2 

DCSAS 

DSSAS 

TTGOV 

TIGOV 

T2GOV 


-  X  2C^/r-  a 

-  2CQ/r  a 


TNDATA 

BODYCM 

ENGNCM 

CUSTCM 

RIDATA 

RTRICM 

R2DATA 

RTR2CH 

QRICM 

QR2CH 

CONTCM 

TRANCM 
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TRANC 


Name:  TRANG(TIM) 

Function: 

calculate  transient  gust  and  control 

General  reference:  section  5*3*^ 

TIM 

time  (dimensionless) 

VELF 

MPSI 

v/n.R 

THDATA 

OMEGA 

RADIUS 

COSPSlO^ 

SINPSl(3d 

0FRTR2 

;i 

reference  rotor 

TRIMCM 

RA1(30} 

RA2(30) 

RTRICM 

RTR2CM 

RWB(3) 

RHT(3) 

RVT(3) 

BODYCM 

RFV(3,3) 

RSF1(3,3; 

RSF2(3,3> 

RHUBlp) 

RHUB2(3) 

1 

) 

MRAl 

ROTATl 

MRA2 

ROTAT2 

RIDATA 

R2DATA 

VGWBV(3^ 

VGHTV(3) 

VGVTV(3) 

vgrtri(3,30,36) 

VGRTR2(3,30,36) 

VQHUBl'' 

VGHUP2(3j 

VPTRAN(5) 

gust  in  wind  axes 

GUSTCM 

OPTRAN 


GCDATA 


0 


A 


rt 


O 


CONTRL 

Name!  CGNTRL(T, PERIOD, C) 

Function:  calculate  transient  control  time  history 
General  reference:  section  5 •3*^ 

Calculates:  C(t)  =  i(l  -  cos2'TTt/T) 

T  time(sec) 

PERIOD  period  T  (sec) 

C  control  C 
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^  a 


o'Cf  ' 


GUSTO 

Name  I  GUSTO(T, PERIOD, G) 

Function:  calculate  uniform  gust  time  history 
General  reference;  section  5,3.^ 

Calculates,  G(t)  =  i(i  -  cos  2 ort/T) 

time  (sec) 

PERIOD  period  T  (sec) 

^  gust  G 
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GUSTO 


Name:  GUSTC(XC,L,L0,G) 

Function:  calculate  convected  gust  wave  stepe 
Genera.l  reference:  section  5 •3*^ 

Calculates:  G(x  )  =  ■Kl  -  cos2'n-(x  -  L-)/L  ) 
g  S  ^ 

XG  distance  y  (ft  or  m) 

g 

L  wavelength  L  (ft  or  m) 

LO  starting  position  (ft  or  m) 

G  gust  G 
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Function:  Performance 

General  reference:  section  5*2.1 


Operating  condition: 

a)  motion:  1st  number  dimensionless,  2nd  number  dimensional 

1)  velocity  =  ft/sec  or  m/sec 

o  o 

2)  dynamic  pressure,  q  =  Ib/ft*^  or  N/m”^ 

3)  weight,  Cy/v-  =  lb  or  N 

4)  body  motion  -  deg/sec,  ft/sec  or  m/sec 

5)  z  =  ft/sec'^  or  m/sec'^ 

6)  =  rpm 

b)  body  orientation  and  controls  in  deg 

Circulation  convergence: 

a)  tolerance,  CG/S  in  C,p/v~  form 

b)  G/E  =  ratio  error  to  tolerance  (-^  1  if  converged) 

Motion  convergence: 

a^  tolerance,  BETA  (etc)  in  deg 

b)  BETA/E  (etc)  =  ratio  error  to  tolerance  (  •^  1  if  converged) 

Airframe  performance:  section  4.2.6 

a)  aerodynamic  loads:  dimensional 

b)  components: 

1)  angles  in  deg 
2^  loads,  q  dimensional 

3)  induced  velocity,  total  velocity  dimensionless 
Gust  velocity:  dimensionless 

System  power:  dimensional  (HP);  number  in  parentheses  is  percent  total  power 

a  ^  nil  mVi  TViuoT  =  V  U 
-  r - 

System  efficiency  parameters: 

a^  gross  weight,  W  =  j.b  or  N  « 

b)  drag-rotor  ~  Dj,  “  (Pj^  +  O/q-rotor  =  D ^yip 

L/D-rotor  =  W/D  ^ 

=)  drag- total  -  dA- total  - 

L/D-total  = 

d)  figure  of  .erit  -  «  -  1  -  P„„„ige^i/Ptotal 


TMDATA 


VEL 

I  TERM 

EmCTN 

ITERC 

EPCIRC 

AFLAP 

0PRTR2 

GRAV 

SIGMA 

RADIUS 

OMEGA 

DENSE 

VELF(3) 

V CLIMB 

VSIDE 

CWS 

HMASS 

NAM 

NDM 

NBMl 

NIMl 

NGMl 

NBM2 

NTM2 

NGM2 

VGWB(3)  gust  in  wind  axes 

VGHT(3) 

VGVT(3) 

VGHUBlO) 

VGHUB2(3) 

vcntrl(ii) 

'^RETFT 

PHIFT 

THETFP 

PSIFP 

THETAT 

PSIT 

DVB0DY(6) 

DOMEGA 

DDZF 

SAVE(3l) 


TRIMCH 


BODYCM 


ENGNCM 

RTRICM 


RTR2CM 


GUSTCM 


oun ion 


QBDCM 


PERF 

LMDAWlpj  WKVICM 

LMDAHlO) 

LHMVlh) 

LHnAW2(3)  WKV2CH 

LHDAH2(3) 

LHDAV2(3) 

BlMS(lO)  CONVGH 

COUNTC 


PERFRl 


Name:  PRRFF1(P,PGFP,PI,PINT,K),PN) 

Function:  calculate  and  print  rotor  performance 
General  reference:  section  5*2.1 


Operating  condition: 


TPP 


1  0  h 

°  1  Ps 

-Pc  -ps  ' 


HP 


o< 


HP 


0 

is 


TPP 


Pc 


^  (^C^CP  ^(^C  ^IsW 
^r^s^CP"  ^(^s  "  ^Ic^HP 


Harmonics  of  gimbal  motion:  section  5*1*2 
Rotor  forces  and  motion: 

shaft  axes  (-S),  tip  path  plane  axes  (-T),  wind  axes  (L  or  X) 
coefficient  (Cx-),  coefficient/solidity  (Cx6-),  dimensional  (x-) 

Rotor  power:  LIIEAL  =  (see  also  section  2.4,3) 


p 

total  power 

PCPP 

climb  and  parasite  power 

PI 

induced  power 

PINT 

interference  power 

PO 

profile  power 

PN 

non-ideal  power 

OPUNIT 

VEL 

MPSI 

MHARM 

MHARHF 

TMDATA 

DENSE 

TRIMCM 

NAM 

BODYCM 

Nm 

ENGNCH 

T75 

CONTCM 

TIC 

TIS 

fz(30,36) 

F_/ac 

AESICM 

ALPHA (30, 36) 

Zf 
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PERFRl 


vind(3,30,36) 

LAMBDA 

(due  to  other  rotor) 

WKVIGM 

vint(3,30,36) 

LAMBDI 

WKV2CM 

RADIUS 

SIGMA 

MRA 

TYPE 

NBLADE 

HINGE 

RIDATA 

MUX 

MUY 

MUZ 

OMEGA 

DRa(30) 

RA(30) 

ALFHP 

PSIHP 

MTIP 

MAT 

N3M 

NIM 

NGH 

NUGC 

NUGS 

RTRIGM 

T75;)LD 

NU(20) 

MDICM 

ETA(2,10) 

WT(ll) 

WTO 

WTC 

WTR 

bending  mode  at  tip 

FHUB(6) 

n'Di  HM 

GLS 

CXS 

BETAO 

BFTAG 

BETAS 

GIRG(36) 

BETA (21, 10) 
THE!.,;  21, 5) 
BETAG(21) 
PHI(10,16) 
PSID(10,6) 

MNRIGM 

QSSTAT(lO) 

PISTAT 

PESTAT 

MNSGM 
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LOAD 


Name :  LOA  D( LEVELl , LEVEL2 ) 

Function:  loads,  vibration,  and  noise 


Airframe  vibration:  section  5*2.8 

Vibration  point  location:  sections  4,1.3,  4.1.5 

LEVELl 

wake  analysis  level  for  rotor  #1 

LEVEL2 

wake  analysis  level  for  rotor  §2 

MHARMF(2) 

TMDATA 

0PRTR2 

FSGALE 

LSCALE 

GRAV 

TRIMCM 

TRATIO 

FSCG 

WLCG 

BLCG 

BDDATA 

NBLDl 

RIDATA 

OMEGA 1 

RTRICM 

NBLD2 

R2D/.TA 

OMEGA 2 

RTR2CM 

MVXRE(3,3) 

MSTAR 

BODYCM 

REn.ER(3,3) 

VELF(3) 

NAM 

THETAT 

PSIT 

CONTCM 

PHI1(10,16) 

HNRICM 

PHI2(10,16) 

MNR2CM 

KVIB 

LA DATA 

FSVIB(IO) 

WLVIB(IO) 

BLVIB(IO) 

ZETA(3,10,10) 
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IjCAItttl 


Name*  LOADRl( LEVEL) 

Function:  calculate  and  print  rotor  loads 


Print  aerodynamics  (function  r  and  H>): 

a)  dimensionless  quantities  generally,  angles  in  dec 
J)  velocity  in  nonrotating  shaft  axes 

c)  interference  induced  velocity  is  that  due  to  other  rotor 
gust  components  in  velocity  axes  other  rotor 


Force/c^^^^  (dimensionless) 


Vc  2^  (c/c  )c  =  L/c 

-  tU  {o/c  )c,  =  d/c 

,  2/  ?  ,'"6^^  d  '  mean 

-  2^  (c  /c  )c  =  m/c 

DR/Q  -  ^  Wn 

UK/u  -  jD  (c/c  )cj  ,  =  D 

Fz/a  -"““r 


D/c 

M/c 


“idUAdi  raaial 
*  '  V‘=„ea„  - 


/o» 


FX/C  =  F  /c 

r  mean 

=  M/c 
a'  mean 


MA/C 
FR/C  =  F^< 
FRT/C  “  F ^c 

Forces  (dimensional) 


'mean 

'mean 


L  »  section  lift 

D  =  section  drag 

M  =  section  pitch  moment 

DR  =  section  radial  drag 

PZ  =  F^  =  dT/dr 
FX  =  F 

V 

HA  =  M 

a 

FR  =  F 

r 

FRT  =  F 

r 

Blade  section  power;  section  5.2.1 
CP/s  =  d(Cp/T-)/dr 
F  =  section  power 


Ib/ft  or  N/m 
Ih/ft  or  N/m 
ft-lb/ft  or  m-N/m 
Ib/ft  or  N/m 
Ib/ft  or  N/m 
Ib/ft  or  N/m 
ft-lb/ft  or  m-N/m 
Ib/ft  or  N/m 
Ib/ft  or  N/m 


HP/ft  or  Hp/m 


LEVEL 

OPUNIT 

HPSI 


level  of  wake  analysis 


TMDATA 


LOAERl 


DENSE 

DPSI 

C0SPSI(36) 

SINPSI(36) 

TRIMCM 

TYPE 

RADIUS 

NBLADE 

OPSTLL 

CH0RD(30) 

inflow(6) 

KRA 

RIDATA 

OMEGA 

CHEAN 

RA(30) 

MUX 

MUY 

MUZ 

NBM 

NIW 

NGM 

PINTER(36} 

pburst(36) 

RTRICM 

ETAT(2,1C) 

ETA(2,10,30) 

bending  mode  at  tip 

bending  mode  ^  i  =  1  to  MRA 

MDICM 

DRV 

WIDATA 

vgust(3,30,36) 

GUSTCM 

GAM(30,36) 

cirg(36) 

QRICM 

HHLOAD 

MALOAD 

MRLOAD 

RL0AD(20) 

NPOLAR 

MWKGMP 

MNOISE 

RANGE(IO) 
ELVATN(lO} 
AZMUTH( 10) 
NPLOT( 75) 

LI BATA 

SAVEM(36,78) 

LOINCM 

M0TI0N(78) 

AEMNCM 
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LOAERl 


3TATE(30,36,3)  AESICH 

DCLMAX?30,36} 

DCDMAX(30,36) 

DCMMAX(30,36j 

MEFFf30,36,3) 

aeff(30,:)6.3) 

dclds(30,36) 

dcdes(30,36) 

DGMDS(30,36) 

save(30,36,19) 

vind(3.30,36)  WKVICM 

LAMBDA 

vwb(3,36} 

vht(3,36) 

vvt(3,36) 

VOFFbi36) 

LAMBDMO) 

LAMBDH(3) 

LAMBDV(3) 

LAMBD0(3) 

V0RH(3,36) 

VINT(3,30,36)  WKV2CM 

LAMBDI 
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LOADHl 


Names  LOADHl 

Functions  calculate  and  print  hub  and  control  loads 


Root  loads s 

MCON  =  Cm 

"'con 

FHUBX  =  Cf^^- 

MHUBX  =  Cmy»- 

FHUBY  =  Cfy<r 

MHUBZ  =  Gm/cr- 

FHUBZ  =  Cf  /j- 

Hub  loads s 

FHUBH  = 

FHUBMX  = 

FHUBY  =  Cy/ct- 

FHUBMY  =  Ch^o- 

FHUBT  =  C^<r- 

FHUBQ  =  CQ/«r- 

1  ^ 

Harmonic  analysis s  f  -  —  ^ 

n  J  j'= » 

F.  K 

J  n 

Dimensional  loads s 

root  force  =  c/r) 

root  moment  =  ^riV(c/R) 

hub  force  =  Ns^n^^Cc/R)  =  (aR)^TTR^r~ 

hub  moment  =  N^r7.^^(c/R)  =  -P  {nR)^'»TR^tr- 


MHARM 

MPSI 

TMDATA 

nbt,ade 

RADIUS 

TYPE 

RIDATA 

CMEAN 

GAMMA 

OMEGA 

NBM 

NTM 

RTRICM 

DENSE 

DPSI 

C0SPSI(36} 

SINPSl(36) 

TRIMCM 

MHARML 

NPL0T(75) 

SENDUR(12) 

CMAT(i2) 

EXMAT(12) 

KFATIG 

for  hub  and  control  loads 

LI DATA 
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LOADHl 


MPAER0(36)  (Mb  /ac 

cmxa(36 
cmza(36 
gfxa(36 
CFZA(36 

cfha(36 

SAVEM(36,78) 

MB 

SB 

10 

SQ(2,10) 

IQA(2,10) 

IFXO 

IHXO 

IP(5) 

IPP(5.5) 

IP0(5) 

IQ0DQ(2,10)  summed  over  q . 

•  J 

SPQ(5.10) 


AEFICM 


LDMNCM 

INCICM 
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"vr 


loads 1 


Name:  LOArsl(R) 

Function:  calculate  and  print  blade  section  loads 
General  reference:  sections  5.2.2,  5.2.3,  5.2.4 

Azimuth  loop:  PHIX  =  T •  t 
PHIZ  = 

T  =8 


FXS-X  =  Gfy<r 

FXS-R  = 

FXS-Z  =  Cf  A- 
■‘z 

CENT  =  Cf 

■^cent' 


HXS-X  = 

=  C^Jrr- 

-  ^'"tors/- 


(-~B  for  shaft  axes,  P  for  principal  axes) 


Harmonic  analysis: 

Dimensional  loads: 

forces 

moments 


=  (*SA)itnVR  -  ^rzV(c/B) 

(!S/a)  =  5nV(o/B) 


radial  station  r/R 

MPSI 

HHARM 

D0FT(4) 


DENSE 

DPSI 

C0SPSI(36) 

SINPSI(36; 

TYPE 

NBLADE 

RADIUS 

MRA 


OMEGA 

CMEAN 

GAMMA 

RA(30) 

DRA(30) 

NBMT 

NBM 

NIM 


TMDATA 

TRTMCM 

RIDATA 

RTRICH 
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LOADSl 


MHARML 

SENDUR(6) 

cmat(6) 

exmat(6) 

KFATIG 

NFlOT(75) 

for  section  loads 

LIDATA 

ETA(2, 10,30) 
DELI 

DEL2 

DEL3 

DEL4 

DEL5 

bending  inodes  at  r^,  i  =  1  to  MRA 

MDICM 

fxaero(30,361 

fzaero(30,36; 

maaero(30,36] 

fraero(30,36; 

1  F  /ac 

1  FVac 

1  M  /ac 

1  F^ac 

AES 1 CM 

BETA(21,10) 

MNRICM 

MB 

• 

• 

« 

IPO 

SAVEM(36,78) 

LDMNCM 
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LOADIl 

Name:  L0ADI1(R,Q,TR,ZR,EPR,ER) 

Function:  calculate  inertia  coefficients  for  section  loads 
General  reference:  sections  5*2.2,  5*2. 3»  5*2.^ 

Blade  pitch:  section  2.3.5 

CS  =  cos©,  SN  =  sine,  TR  =  e(r) 

W  =  (z?-xk),  WP=(z^-xk),  WPP  =  (zT-xk) 

'o  o  _j^  o  o'  'o  o 

WXI  =  (zT!-x‘^  -  x^-k ) 
o  o  I 

ZR  =  li(r),  ER  =  ^.(r),  EPR  =  ^/(r) 

WR-(zlt-xk),  ,  ,  WPR  =  (z^-xk),  ,  ,atr 
^  o  o  -^trim  o  o  '^trim’ 

WRXC  =  (z'^-xTc  -  x  Jc) ,  at  r 
'  o  o  Q' 

EPXIO(NBM)  =  (  ^^-k  Xj)  at  r  =  e 
CE(NBM)  --  x^Tc-  Xj^)  d^ 

CMR(MRM+l)  =  -  r)  m  d 

UFA  =  (z  I*-  X  k),  WPFA  =(z^-xk)  at 

'0  0  0  0  rA 

X  =  \(  ^).  XR  =  if^(r) 


ZR(5) 

EPR(2,10) 

ER(2,10) 

DEBUG 

T75 

EFLAP 

ELAG 

XFA 

RFA 

ZFA 

RCPL 

NOPB 

MRM 


radial  station  r/R 

mean  deflection  q . 

0 

pitch  0^  at  r 
Ik  at  r 
Ik  at  r 
at  r 


TMDATA 

CONTCM 

RIDATA 
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LOADIl 


NBM  RTRICM 

Nm 

NGM 

NBMT 

MASS(51) 

ITHETA(51) 

XI(51) 

TWIST(51) 

ETA(2,10,51)  bending  modes  at  r=(j-l)Ar,  j=l  to  MRH+1  MDICM 

ETAP(2,10,51) 

ETAPP(2,10,5l) 

ZETA(5,51)  torsion  modes  at  r=(j-l)2vr,  j=l  to  MRH+1 

ETAPH(2,10) 

EFa(2,10)  bending  modes  at  r  =  r„ 

EFAP(2,10) 

DELI 

DEL2 

DEL3 

DEL4 

DEL5 

MB  LDMNCH 


LOADF 


Namej  LOADF(S,MPSI,K, SB, C.H,  DAMAGE, SEQ) 
Function:  calculate  fatigue  damage 
General  reference:  section  5*2.9 


Input: 

S(MPSI) 

HPSI 

K 

SE 

H 

C 


vector  of  load  S.,  j  =  1  to  MPSI;  dimensional 
number  of  azimuthal  stations;  maximum  36 
parameter  K  in  fatigue  damage  calculation 
endurance  limit  Sg  (dimensional) 
material  exponent 
material  constant 

Q 

S-N  curve  approximated  by  N  =  - — 

(S/Sg  -  1)" 


Output: 

DAMAGE  damage  fraction  per  rev  (only  calculated  if  S_  >  0, 

C  >  0,  and  M  O) 

SEQ  equivalent  i  peak-to-peak  load  (only  calculated  if 

MM) 
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LOADM 


Name:  L0ADM(F,MPSI,FME.AN.FHPP) 

Function:  calculate  mean  and  half  peak-to-peak 

Input: 

F(MPSI) 

MPSI 

Output: 

FMEAN  mean  load 

FHPP  i  peak-to-peak  load 


load  F . ,  j  =  1  to  MPSI 
«3 

number  of  azimuthal  stations 
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D 


GEOHPl 


LEVEL 

waJte  analysis;  for  prescribed  wake, 

2  for  free  wake  geometry 

MPSI 

TYPE 

TMDATA 

RIDATA 

MWKGMP 

JWKGMP(8} 

NWKGMP(4) 

LI DATA 

KFW 

KDW 

KNW 

KRW 

KRWG 

WIDATA 

KFWG 

GIDATA 
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PCLRPP 


Name :  POLRPP( A ,  f®  A , RA , MPSI , ISUB , NPLOT , DA , NUPP) 
Function:  printer-plot  of  polar  plot 


(J-NCoL*  fol  +  KC.oL_  = 


A 

MRA 

ra(mra) 

MPSI 

ISUB 

NPLOT 

DA 

NUPP 


array  to  be  plotted 

number  of  radial  stations 

radial  stations  r^^,  i  =  1  to  MRA 

number  of  azimuthal  stations  't.  =  j  , 
j  =  1  to  MPSI,  360/HPSI  ^ 

first  dimension  of  array  A;  positive  if  first  subscript 

^  ^  **S^cx  i/xvc  XX  axx5u  5  uuc?  Cx  x  u  u  XO 

n;  data  plotted  every  n-th  step 

plot  increment:  last  digit  of  integer  part  of 
A/DA  is  plotted  (if  multiple  of  NPLOT) 

unit  number  for  printed  output 
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HISTPP 


Name  .*  HISTPP(  A ,  MRA ,  RA ,  MPSI ,  ISUB ,  NPLOT ,  NAME ,  NUPP) 
Function:  printer-plot  of  azimuthal  time  history 


let  c  minimum,  d  -  maximum  values  over  azimuth 
l)  d>0,  c<-.03d  or  c^O,  d'>.031c| 


use  b  = 

[^max(d,  lc\)J 

\  A 

-k  ' 

2) 

d  >  2jcl  , 

r  c  >  -.03d 

use  b  = 

CO 

1  A 

0 

3) 

c  <  -  2  |d| 

1  ,  d  <  .03  c 

use  b  = 

[uo 

-W 

4) 

otherwise. 

use  mean  =  [i(c  +  d)] 

and  b  - 

[  max  (  mean-c  ,  d-mean  )  J 

1  M 

v«.eo.y\-.V5 

mean  =  AH  =  KM  * 10**NM 
b  =  B  =  K  *  10^‘*N 


to  convert  F  to  K*10^  (K  =  1  to  9) 

a)  if  F  =  0,  then  F  =  .99 

b)  N  =  [log  |Fl] 

if  F  <  1.,  then  N=  N  -  1 

c)  K  =  [iF!  /10**n1  +  1 

if  K  =  10,  then  N  =  N  +  1  and  K  =  1 
if  F  <  0,  then  K=  -K 

d)  F  =  K  » 10**N 


HISTPP 


A 

MRA 

ra(mra) 

MFSI 
I  SUB 
NPLCiT 

NAME 

NUPP 


array  to  be  plotted 

secondary  variable:  number  of  values  (minimum  l) 

secondary  variable:  values  r. ,  i  =  1  to  MRA; 
alphanumeric  labels  if  NPLOT  LT  0;  not  used  if 
MRA  EQ  1 

number  of  azimuthal  stations  «  j  =  1  to  MPSI, 

=  360/MPSI  ^ 

first  dimension  of  array  A;  positive  if  first  subscript 
is  r^,  negative  if  first  subscript  is 

number  of  values  of  secondary  variable  per  plot; 
minimum  1  and  maximum  3s  negative  for  alphanumeric 
labels;  not  used  if  MRA  EQ  1 

name  of  secondary  variable,  4  characters;  not  used 
if  MRA  EQ  1 

unit  number  for  printed  output 
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NOISRl 


Name:  N0ISR1(RANGE,ELVATN,AZMUTH) 

Function:  calculate  and  print  far  field  rotational  noise 
General  reference:  section  5 ‘2. 10 

Calculate  constants:  GSTR  =  cos^^(l  -  M^) 

FT  =  -N' 

FD  =  N^/4ir^(l  -  M^) 

FL  =  -N^jnsin9^/4TrCgC^(l  - 

FR  =  -N^flcos©  ATrc^«7;;(l  -  M 
IT  s  3r 

FB  =  Nflcos©  /c  (1  -  M  ) 
r  S  DC 

FS  =  Nllg^/Cg 


Harmonic  analys 

1  -J  ^4^1 

is  of  loads:  F  =  ^  F.e  J  K 

n  J  jsi  0  n 

RANGE 

range  s  (dimensional) 

ELVATN 

elevation  €5  (deg) 

AZMUTH 

azimuth  (deg) 

MPSI 

OPUNIT 

TMDATA 

DPSI 

DENSE 

CSOUND 

cospsi(36} 

SINPSI(36) 

TRIMCM 

OMEGA 

CMEAN 

MUX 

MUY 

MUZ 

RA(30) 

DRA(30) 

RTRICM 

NBLADE 

CH0RD(30) 

SIGMA 

RADIUS 

MRA 

TYPE 

RIDATA 

AXS(30j 

opnoisAJ 

MHARMN(3) 

MTIMEN(3) 

A  /c^ 
xs' 

LI DATA 

-l6l- 


NCISRl 


FXA(30,36)  F  /ac 

FZA(30,36)  FVac 

FRA(30,36)  F^ac 

BETAC 

BETAS 


AESICM 


QPICM 
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BESSEL 


Name:  BBSSEL(NB,XB,BJ) 

Function:  calculate  J  Bessel  function 

Input: 

NB  order  of  Bessel  function,  n 

XB  argument  of  Bessel  function,  x 

Output: 

BJ  Bessel  function  J  (x) 
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RAMF 


Name:  RAMF(LEVELl,LEVEa:,2,0PLMDA) 

Function:  calculate  rotor/airframe  periodic  motion  and  forces 
General  reference:  section  5=1 ‘13 


Test  motion  convergence:  section  ‘^.1.4 
Test  circulation  convergence:  section  5*1 ‘12 


LEVEL 1 
LEVEL2 


integer  parameter  specifying  rotor  and  rotor  it2 
wake  analysis:  0  for  uniform  inflow,  1  or  2  xoT 
nonuniform  inflow 


CFLMDA 


integer  parameter:  0  to  suppress  inflow  update 


MPSI 

MHARH(2) 

MHARMF(2) 

I  TERM 

EPMCTN 

ITERC 

EPCIRC 

DEBUG 

MREV 

MPSIR 

0FRTR2 

NAM 

NDM 


THDATA 


TRIHCM 
BODY  CM 
ENGNGM 


CMEANl 

NBMl 

NTMl 

NGMl 

CMEAN2 

NBM2 

NTM? 

NGM2 

Bl(21,10) 

Tl(21,5) 

BGl(2l) 

P1(10,16) 

PS1(10,6) 

B2(21,10) 

T2(21,5) 

BG2(21) 

P2(10,16) 

PS2(iO,6) 


RTRICM 


RTR2CM 


MNRICM 


MNR2CM 
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BlMS(lO) 


CONVCM 


COUNTC 

CIRC1(36) 

GTl 

CMXl 

CMYl 

CIRC2(36) 

CT2 

CMX2 

CMY2 

SIGMAl 
SIGMA 2 


QRICM 


QR2CM 


RimTA 

R2DATA 


MODEl 


Name :  MODEl 

Function :  blade  modes 

T750LD 

MDICM 

NBMOLD 

NTMOLD 

DEBUG 

TODATA 

HINGE 

RIDATA 

EFMOIE 

NBM 

BTRICM 

NTM 

T75 

CONTCM 
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MOEECl 


Names  MO DEG 1 

Functions  initialize  blade  mode  parameters 

Linearly  interpolate  data  for  bending  mode  calculations  section  2.3 •! 

Tip  masss  section  2.2.19 

Evaluate  centrifugal  force  for  bending  mode  calculations  section  2.3.1 

CENT  =  di^ 

Linearly  interpolate  data  for  torsion  mode  calculations  section  2.3.3 
Evaluate  pitch  inertia  and  control  system  stiffnesss  sections  2.2.9»5»1*3 

MRS  RIDATA 

MTIP 
XITIP 
EFLAP 
ELAG 
RFA 
RADIUS 
MRH 
FTO 
FTC 
FTR 
WTIN 
VTIPN 
KTOI 
KTCI 
KTRI 
HRI 

RI(51 
XI(51 
XG(51) 

KP2(5l) 

MASsbl) 

ITHETA(51) 

Em(5i) 

Eizz(5l) 

TWIST^l) 

DEBUG  TMDATA 


-167- 


MOIECl 

13  RTRICM 

OMEGA 

EIXXB(5l) 

EIZZBfsij 

HASSB(51) 

TWISTB^l) 

CENT(5l) 

ITHETB(51) 

GJB(5l) 

MASSl(5l) 

ITHETl(5l) 

XII(51) 

XCI(51) 

TWISTl(5l) 

KP2I(51) 

I PITCH 
KTO 
KTC 
KTR 


HOPES i 


Name:  MODEBl 

Function:  calculate  blade  bending  modes 
General  reference:  section  2.3*1 

4  ^11  f  'I  , 

V  •  W  •‘'- 

Co  ''K'  Sl  -wAr 

(■  CEX/a*fc-T’C-  j-r 

X 

Normalize  eigenvector  solution:  using  Galerkin  modes  from  last  call, 

which  was  at  r  =  1 


Blaide  pitch:  section  2.3*5 
Calculate : 

DS  = 

C  = 
DC  = 

B  = 

A  = 


T75 

CONTCM 

DEBUG 

TMDATA 

NOPB 

RIDATA 

RCPL 

KFLAP 

KLAG 

EFLAP 

ELAG 

RADIUS 

RCPLS 

TSPRNG 

RFA 

RPB 

NCULb 

MRB 

NONROT 

HINGE 

MRA 

RROOT 

MRM 

NU(20) 

MDICM 

NUNR(20) 

ETA(2,10,96) 

ETAP(2,10,96) 

ETAPP(2,10,96) 

ETAPH(2,10) 
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MO  DEB  1 


MASSfSl)  inertial  and  structural  data  at  RTRICM 

EIXX(5lj  r  =  e  +  (j-l)Ar,  r  =  1  to  MRB  +  1 

BIZZ(5l) 

TWIST^l) 

CENT(51) 

OMEGA 

NBM 

RA(30) 
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MO  DBG 

Name i  HOnEG(R,EFLAP. BLAG, NCOLB, HINGE, F, 

IF, DDF) 

Function: 

calculate  Galerkin  functions 

for  bending  modes 

General  reference:  section  2.3.1 

R 

radial  station  r/R 

EFLAP 

flap  hinge  offset  e^/R 

BLAG 

lag  hinge  offset  e^/R 

NCOLB 

number  of  functions 

HINGE 

integer  parameter:  0  for  hinged  blade,  1  for 
cantilever  blade 

F( NCOLB) 

Galerkin  functions  f^^ 

DF( NCOLB) 

Galerkin  functions  f^ 

DDF  (NCOLB) 

Galerkin  functions  f" 
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MODEM 


Najne;  MOIEAi 

Function:  calculate  arxiculated  blade  flap  and  lag  modes 
General  reference:  section  2.3 i 2 

f*  2 

Calculate:  F  =  '^v^mdr,  G  =  jvjmdr 

DEBUG 

MRB 

EFLAP 

ELAG 

KFLAP 

KLAG 

RADIUS 

MRM 

RFA 

RPB 

MRA 

RROOT 

RA(30) 

OMEGA 

NBM 

MASS(5l)  section  mass  at  r=e  +  (j-l)Ar,  j  =  1  to  HRB+1 

NU(2C) 

NUNR(20) 

ETA(2,10,96) 

ETAP(2,10,96) 

ETAPP(2,10,96) 

ETAPH(2,10) 


TMDATA 

RIDATA 


RTRICM 


HDICM 


MODETl 


Name:  MODETl 

Function:  calculate  blade  torsion  modes 
General  reference:  section  2.3.3 


Bvalvate  Galerkin  functions  at  r: 
Calci late : 

A  = 

B  = 

C  = 


X  =  ix{r  -  rp^)/(l  -  rp^) 

CC.-3-/si^g.*-y 

‘'O 


w  v: 


Normalize  eigenvector  solution:  using  Galerkin  functions  from  last 

iteration,  which  was  at  r  =  1 


DEBUG 

MRB 

RFA 

RADIUS 

MRM 

NCOLT 

MRA 

IPITCH 

KTO 

KTC 

KTR 

OMETA 

NTH 

RA(30) 

ITHETA(5l)  le  at  r  =  r„  -Ko-OAr,  j  =  1  to  MRB+1 

GJ(5l)  GJ  at  r  =  rpf+(j-l)Ar;  j  =  i  to  HRB+l 

WT(ll) 

WTO 

WTC 

WTR 

ZBTA(5,92) 

ZETAP(5,92) 


TMDATA 

RIDATA 


RTRICM 


MDICM 
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MOEEKl 


Name :  MODEKl 

Function:  calculate  kinematic  pitch-bending  coupling 
General  reference:  section  2.3.4 


DEBUG 

TMDATA 

T75 

CONTCM 

PHIPL 

RIDATA 

.HIPH 

RPH 

RPB 

XPH 

KPIN 

DEL3G 

ATANT{P(10) 

ETA(2,10) 

benc’ing  modes  at  rpg 

MDICM 

ETAP(2.10) 

KPB(lO) 

KPG 

NBM 

RTRICM 
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MODEDl 


Name:  MODEDl 

Function:  calculate  blade  root  geometry 
General  reference:  section  2.2.1 


DEBUG 

T75 

CONE 

DROOP 

SWEEP 

FDROOP 

FSWEEP 

DELI 

DEL2 

DEL3 

DEL4 

DEL5 


IMDATA 

CONTCM 

RIDATA 


MDICM 
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INRTCl 


la 


Name:  INRTCl 

Function:  calculate  blade  inertia  coefficients 
General  reference:  section  2.2.19 


Blade  pitch: 
Calculate: 


section  2.3*5 

CS(MRM-i-l)  =  cos©-,  SN(MRM+1)  =  sin© 
CM(HRM+l)  = 

CMR(MRM-t-l)  =  Vf 

CMRR(MRM+l)  =  Cf 

CXIM(MRM-i'l)  =  ‘)r  Xx-^‘5>e> 

CXIPM(MRM+l)  =  V  0 

DEM(NBM,MRM+l)  =  Vr 

DERM(NBK,MRM-^1)  = 

CEPEP(NBM,NBMT,MRMfl)  =  C  1'^  '  f  i 
X(2,NTM,NBMT,MRM-<-l)  =  X 


kj 


a)  X  = 

b)  XH  = 

c)  X  = 
XCFA  =  at 

XCE  =  x^  at  e 
XIE  =  Xj  at  e 
KP2'IVP  =  kp 


u  ^ 

X,  .  for  k  1  and  k  =  0 
kj 


DEBUG 

T75 

MRM 

NOPB 

RCPL 

RFi 

ZFA 

XFA 

ELAG 


TMDATA 

CONTCM 

RIDATA 
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INRTCl 


RADIUS 

MBLAIE 

MRA 

EFLAP 

IB 

NBH 

NTM 

NGM 

NBMT 

RA(30) 

I PITCH 
MASS(5l) 
ITHETA(51) 
XI(51) 

XC(51) 

KP2(5l) 

TWIST(5l) 

ETA(2,10,51) 

ETAP(2,10,5l) 

ETAPP(2,10,5l) 

ZETA(5,51) 

ZETAP(5,51) 

EFA(2,10) 

EFAP(2,10) 

ETAPH(2,10) 

DELI 

DEL2 

DEL3 

DEL4 

DEL3 

MB 


inertial  data  at  r=  j  =  1  to  MRM+1 

bending  modes  at  r  =  j  =  1  to  MRM+1 
torsion  modes  at  r  =  (j-l)Ar,  j  =  1  to  MRM+1 
bending  modes  at  r^^ 


RIDATA 


P.TPICM 


INC 1 CM 


XAPQ(2,5,4,30) 
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HOIEPl 


Name :  MODEPl 

Function;  print  blade  modes 


TYPE 

HINGE 

NCOLB 

NONROT 

NCOLT 

RCPL 

EFLAP 

ELAG 

KFLAP 

KLAG 

RCPLS 

TSPRNG 

RADIUS 

Ot-IEGA 

NBM 

NTM 

NGM 

NUGC 

NUGS 

KTO 

KTC 

KTR 

IB 


bending  modes  at  r  =  ( j  =  1  to  11 


torsion  modes  at  r  =  (o-l).l»  j  =  1  to  11 


MB 

SB 

IP(5) 

T750LD 
NU( 20) 

NUNR(20) 

ETA(2,10,11) 

ETAP(2,10,11) 

ETAPP(2,10,ll) 

WT(ll) 

WTO 

WTC 

WTR 

ZETA(5,11) 

2ETAP?5.1l) 

ETAPH(2,10) 

KPB(IO) 

KPG 


RIDATA 


RTRICM 


INCICK 


MDICM 


MODEPl 


DELI  MDICM 

DEL2 

DEL3 

DEL4 

DBL5 
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I 


BODYC 


Name :  BODYC 

Function:  initialize  airframe  parameters  at  trim 

Wind  tunnel  trim  case:  section  4.1.3 

r,  Rgp  with  rotations:  sections  4.1.3,  4.1.5 

Free  flight  trim  case:  section  4.1.1 
Calculate  R  :  section  4.2.1 

Calculate  R^  I*Rg,  -M*(V  x  )Rg,  G,  (Vx)Rgkp:  section  4.2. 
Airframe  gust  velocity  in  body  axes:  section  4.1.4 


THETFT 

PHIFT 

PSIFP 

THETFP 

THETA T 

PSIT 

DEBUG 

VEL 

OPTRIM 

MSTAR 

MSTARG 

ISTAR(3,3) 

RSF10(3,3) 

RSF20(3,3) 

RHUB10(3) 

RHUB20(3) 

RWB0(3j 

RHT0(3) 

RVT0(3) 

R0FF0(3) 

RSF1(3,3) 

c.\J  j 

RHUBl(3) 

RHUB2(3) 

RWB(3) 

RHT(3j 

rvtO) 

R0FF(3) 

VXRBKF(3) 

MVXRE(3,3) 

GMTRX(3,3) 

IBODY(3.3) 

REULBRb.3) 

RFV(3,3) 


CONTCM 

TMDATA 

BODYCM 
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BODYC 


HFB(3,3)  BODYCH 

KE(3) 

VELF(3) 

V CLIMB 
VSIDB 

VGWBVf3}  CUSTOM 

VGHTV(3) 

VGVTV(3) 

vgwbf(3J 

VGHTF(3) 

VGVTF(3) 
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ENGNC 


Name:  ENGNC 

Fvinction:  initialize  drive  train  parameters  at  trim 

Engine  damping:  section  4.3 •! 

Drive  system  inertia:  section  5*3 

Drive  system  spring,  damping,  mass  matrices:  section  5* 1*9 
Drive  system  static  elastic  matrix:  section  5«1*10 
Calculate  C^  :  section  5.1 *5 
Calculate  C^  :  section  5*1*9 


DEBUG 

OPENGN 

TODATA 

0PRTR2 

TRIHCH 

NBLDl 

RIDATA 

NBLD2 

R2DATA 

IBl 

OMEGA 1 

gamma 1 

CD1(2) 

CPSI1(2) 

RTRICM 

IB2 

0MEGA2 

GAMMA 2 

CD2(2) 

CPSI2(2) 

RTR2CM 

101 

QTl 

QDZl 

I NCI CM 

102 

QT2 

QDZ2 

INC2CM 

CQSl 

QRICM 

nnoo 

^  -if  on  3 

KIGOVE 

KIGOVl 

KIG0V2 

GSE 

GSI 

KEDAMP 

ENDATA 
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ENGNC 


QTHRTL  ENGNCM 

lENG 

IMIl 

KMI2 

KMR 

HKEl 

KHE2 

KPGOVE 

KPGOVl 

KPGOV2 

TIGOVE 

TlGOVl 

T1G0V2 

T2G0VE 

T2GOV1 

T2GOV2 

QEiWMP 

IRSTAR 

MENG(6,6j 

seng(6,6) 

DENG(6,6) 

HENG0(2,2) 
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MOTNCl 


Name:  MOTNCl 

Function:  initialize  rotor  parameters  at  trim 

Calculate  ^at'  sections  2.4.2,  4.1.2 

Calculate  section  4.1.4 

Rotor  gust  velocity  in  shaft  axes:  section  4,1,4 
Calculate  c,  c:  section  4.2.2 
T 

Calculate  c  «  section  4.2.5 
Calculate  section  4.1.2 


DEBliG 

MPSI 

NSCALE 

ISCALE 

FSCALE 

LSCALE 

IB 

OMEGA 

MTIP 

MUX 

MUY 

MUZ 

ALFHP 

PSIHP 

HAT 

RGUST(3,3) 

CHUB(6,16) 

CBHUB(3,3) 

CHUBT(16,6) 

ROTATE 

NBLADE 

RADIUS 

MRA 

NEM 

DVL0DY(6) 

VGUSTV(3,30,36)  gust  at  rotor  disk,  velocity  axes 
VGUSTS(3 ,30,36)  gust  at  rotor  disk,  shaft  axes 
VGUSTH(3)  gust  at  rotor  hub,  velocity  axes 


THDATA 

TRIHCM 


RTRICH 


HI DATA 


BDDATA 

CONTCM 

GUSTCM 
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HOTNCl 


VELF(3) 

RFV(3,3) 

REULER(3,3) 

RSF(3,3) 

RHUB(3) 

AM0DE(6,10) 


BODYCM 
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BODYMl 


Name:  BODYMl 

Function: 

calculate  airframe  transfer  function  matrix 

General  reference:  section  5. 1.8 

DEBUG 

TMDATA 

DOF(l6) 

MHARMP 

airframe  degrees  of  freedom 

FSCALE 

TRIMCM 

NBLADE 

RIDATA 

OMEGA 

RTRICM 

DPSI21 

(rad):  O.  for  rotor  jj^l 

chubt(i6,6) 

AMASS(lO) 

BODYCM 

ADAMPS(10; 

ASPRNG(10 

ADAMPA(10 

IBODY(3.3 

mvxre(3,3 

GM'mX(3,3; 

MSTAR 

NAM 

HB0DY(16,6,10) 

RHICM 
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ENGNMi 


Name :  ENGNMI 

Function:  calculate  drive  train  transfer  function  matrix 
General  reference:  section  5«1*9 


DEBUG 

MH^RMF 

IX)F(6) 

drive  train  degrees  of  freedom 

TMDATA 

FSCALE 

TRIMCM 

NBLADE 

RIDATA 

OMEGA 

DPSI21 

CD(2) 

(rad);  0.  for  rotor 

RTRICM 

MENG(6,6) 

SENG(6,65 

DENG(6,6) 

Nm 

ENGNCM 

HENG(6,10) 

RHICM 
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WAKEUj 


Name:  WAKEUl 

Function:  calculate  uniform  wake-induced  velocity 
General  reference:  section  2.4.3 

Lagged  thrust  and  moment:  section  5.1.12 

Vectors  for  aerodynamic  interference:  section  4.2.6 

Interference  induced  velocity:  section  4,2.6 

DEBUG 

OPGRND 

HAGL 

MPSI 

DPGX 

C0SPSI(36} 

SINPSI(36) 

LSCALE 

FSCALE 

MRA 

RADIUS 

ROTATE 

FACTOR 

KHLMDA 

KFLHDA 

FXLMDA 

FYLMDA 

FMLMDA 

KINTH 

KINTF 

Kl^TWB 

KINTHT 

KINTVT 

I^^FL0W(6) 

RA(30) 

OMEGA 

MUX 

KUY 

MUZ 

MRAO 

RADUSO 

OMEGAO 

R3F(3,3) 

RHUB(3) 

RWB(3) 

RHT(3) 

rvtO) 

KE(3) 


TMDATA 


TRIMCM 


RIDATA 


RTRICM 


R2DATA 

RTR2CM 

BODYCM 


WAKEUl 


CT 

Cip 

QPICM 

CMY 

CMv 

CMX 

CTOLD 

nMV/-\T  TV 

WKVICM 

CMTOLD 

vind(3,30.36) 

LAMBDA 

FGS 

COSE 

ZAGL 

vint(3,30.36) 

LAMBDI 

LAMBDtf(3) 

lambdhO) 

LAMBDV(3) 

LAMBD0(3) 

EINTW(3} 

EINTH(3) 

e:ntv(3) 
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WATilNl 


Name:  WAKSNl (LEVEL) 

Function:  calculate  non-uniform  wake  induced  velocity 
General  reference:  section 


Calculate  section  3*1  *3 

^TF  ^TS^SF 


^21  "  ^^SF^other  rotor^TF 


Lagged  circulation:  section  5 •1*12 

Interpolate  induced  velocity:  linear  interpolation  between  inflow 

points,  constant  beyond  first  or  last 
point 

Calculate  mean  induced  velocity:  TPP  normal  component,  area-weighted  mean 


LEVEL  rotor  wake  level:  0  for  uniform  inflow  (only 

replace  old  circulation) 


DEBUG  TtiDATA 

MPSI 

DPSI  TRIMCM 

MRA  RIDATA 

ROTATE 

INFL0W(6) 

RA(30)  RTRICH 

draOo) 

DP21M  (rad)j  0.  for  rotor  #1 

DPSI21  (rad);  for  rr  tor  #2 

MRAO  other  rotor  R2DATA 

ROTATO 

RA0(30)  RTR2CM 

DRAG (30) 

NG(30)  W2DATA 


MRG 

NL(30) 

MRL 

FACTOR 

OPVXVY 

KNW 

OPRTS 

NLO(30)  other  rotor  W2DATA 

MRLO 

RSF(3,3)  30DYCM 

RSFO ( 3  .3 )  0  ther  ro  tor 
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WAKENl 


GAHf30,36)  QRIGM 

CRC(36) 

BETAC 

BETAS 

BETACO  other  rotor  QR2GM 

BETASO 

GAM0LD(30,36)  WKVICM 

CRC0LD(36) 

vind(3,30.36) 

LAl-ffiDA 

vintO,30,36) 

vorh(3,36) 

LAMBDI 

vwb(3,36} 

vht(3,36) 

vvt(3,36) 

VOFF(3,36) 

LAMBEW(3) 

LAMB Dm  3) 

.jAMBDV(3) 

LAMBDO(3) 

MR  WKClCM 

ML 

MI 

MW 

MH 

MV 

MO 

C(3, 20000) 

CNW(3,20000) 


-191- 


INRTOl 


Name:  INRTHl 

Function:  calculate  rotor  transfer  function  matrix 
General  reference:  section  5*1 ‘6 


Aerodynamic  spring  and  damping:  section  2.2.20 

DEBUG  TMDATA 

DOF (15)  rotor  bending  and  torsion  degrees  of  freedom 

D0FT(4) 

MPSI 

HHARH 

RA(30)  RTRICM 

draOo) 

CMEAN 

HUZ 

NUGC 

NUGS 

CGC 

CGS 

GLAG 

CTO 

CTC 

CTR 

NBM 

NTH 

NGM 

NBHT 

GAMA  ^ 


KEPSI(21,36)  TRIHCM 

HRTR(16,16,21)  RHICM 

CT  0;?1CM 

LAMBDA  "  WKVICH 

BETA(21,10}  HNRICM 

THETA(21,5) 

BETAG(21) 

F0RCE( 16,36)  AEFICH 

NBLAEE  RIDATA 


GSB(lO) 

GST(5) 

MRA 

CHORD(30) 

SIGMA 

XA(30) 

XAC(30) 
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INTRMl 


mj(20) 

BTAPH(2,10) 

KPG 

KFB(IO) 

AETA(2,10.3) 

AZETA(5,30) 

WT(ll) 

WTO 

WTC 

wra 

HE 


bending  modes  at  r  ,  i  =  i  to  mrA 
U)rslon  modes  at  i  =  i  to  HRA 


XARi(2,5,4.30) 

KCiDQ(lO.lO) 

• 

MPP(5,5) 

IQDQS(10,10) 

• 

SP(IS(5,10) 


MDICM 


INCICM 
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INRTI 


Name:  INRTl(HX,H,KEEP,LHINV,MMIHV) 

Function:  calculate  inverse  of  transfer  function  matrix 


MX 

H(MX*HX) 

KEEP(MX) 

LMINV(MX+l) 

mminv(mx+i) 


dimension  of  H 

n 

complex  matrix  to  be  inverted 

integer  vector  designating  degrees  of  freedom 
to  be  retained;  0  for  unused  degrees  of  freedom 

scratch  vector 

t cratch  vector 


-194- 


MOTNHl 


Name :  MOTNHl 

Function:  calculate  harmonics  of  hub  motion 

General  reference 

:  sections  5*lt5f  5*1*11 

DEBUG 

TMDATA 

MHARM 

HHARMF 

GRAY 

FSCALE 

LSGALE 

TRIMCM 

RADIUS 

ROTATE 

NBLADE 

RIEATA 

OPHVIB(3) 

OMEGA 

CHUB (6, 16) 

RTRICM 

cbhub(3,3) 

CPSI(2) 

^*f'y(rad)5  0,  for  rotorffl 

DPSI21 

KMASTG(IO) 

KMASTS(IO) 

BODYCH 

RSF(3,3) 

KE(3) 

NAM 

NDfi 

ENGNCM 

DVB0DY(6) 

DOMBGA 

CONTCH 

QSSTAT(IO) 

PISTAl 

MNSCM 

PHI(10,16) 

MNRICH 

PSID(10,2) 

THTG(lO) 

f  vi->  u>_\ 

V  ^ 

PHIO  ('10,16) 
PSID6(10,2) 
THTGO(IO) 

^  (due  to  other  rotor) 

MNR2CM 

ALF(10,6) 

• 

• 

MNHICM 

« 

DPSISO 
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HOTNRl 


Name:  HOTNRl (JSTART) 

Function:  calculate  harmonics  of  rotor  motion 
General  reference:  sections  5 A. 6,  5* 1*13 

Lag  damper  moment:  section  2.2. iS 

Calculate  coning  and  tip-path  plane  tilt:  section  3.1.3 
Calculate  hub  reactions:  section  5.1.7 


JSTART  azimuth  index 

MPSI 

MPSIR 

DEBUG 

MHARH 

MHARHF 

D0FT(4) 

JfBLADE 

GAMMA 

NBH 

NTH 

NGM 

NBMT 

GLAG 

HLD 

DZLD 

CGC 

CGS 

NUGS 

NUGC 

KFB(lO) 

KPG 

ETAPH(2,10) 

ETATIP(2,10)  bending  mode  at  r  =  1 

BO 

BC 

BS 

BETA(2i  ,10) 

THETA(21,5) 

BETAG(21) 

DPSI 

C0SPSI(36) 

SINPSI(36) 

KEFSI(21,36) 

HRTR(16,16,21) 


TMDATA 


RIDATA 

HTRICM 


MDICM 

QRIGH 

MNRICM 

TRIMCM 

RHICM 
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MOTNBl 


Namei  MOTNBl (PSl) 

Function:  calculate  blade  and  hub  motion 

General  reference:  section  5«1«5 

Rigid  pitch  p^:  section  5*1 *3 

PSI 

q(io) 

• 

• 

AEMNCH 

• 

DTT 

HHARM 

TMDATA 

HHARFM 

NBLADE 

RIDATA 

NBM 

RTRICH 

NTM 

NGM 

KFB(lO) 

MDICM 

KPG 

T75 

CONTCM 

TIC 

TIS 

BETA(21,10} 

HNRICH 

THETA(21,5) 

BETAG(2l) 

ALF(10,6) 

• 

MNHICM 

• 

DPSISO 
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AEROFl 


Name :  AEROFl ( JPSI , QT , MQ , HP , CMX , CMZ , CFX , CFZ , CFR ) 

Function:  calculate  blade  aerodynamic  forces 

Calculate  XAP  =  Xy^j^:  section  2.2.19 

Section  velocity  components:  section  2.4.2 
Calculate  U,  M,  <)>  ,  :  section  2.4.1 

4>  in  rad,  in  deg 

Calculate  <vc/V:  section  2.4.7 

Calculate  cos_A_ :  section  2.4.6 

REVFLW  =  1  if  just  crossed  reverse  flow  boundary 

Tip  loss  correction:  section  2.4.5 

Section  forces  and  pitch  moment:  section  2.4.1 

FZ  =  F  /ac  ,  FX  =  F  /ac  ,  FR  =  F  /ac  ,  MA  =  M  /ac 
z'  m’  x  m  r  m  a'  m 

Circulation:  section  2.4.9 

Unsteady  lift,  moment,  and  circulation:  sections  2.4.8,  2.4.9 

LUS  =  L  /ac,  HUS  =  M  /ac,  GUS  =  P  /ac 
us'  *  us'  '  us' 

Maximum  circulation  outboard  r^  :  section 

'-•max 


JPSI 

azimuth  index  j 

QT(4) 

^jtrim 

MQ(IO) 

Mqj^aeroA® 

MP(5) 

^Pj^aero/®-® 

CMX 

CMZ 

Cmg/'^a 

CFX 

Cfj/«-a 

CFZ 

CFR 

Q(10) 

DQ(IO) 

DDQ(IO) 

P(5) 

DP(5)^ 

DDP(5) 

BG 

DBG 

DDBG 

ahub(6) 

dahub(6) 

DDAHUB(6) 

Gf^v-a 

AEMNCM 
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AEROFl 


PS 

DPS 

DDPS 

aemngm 

DEBUG 

MPSI 

TODATA 

DPSI 

FSGALE 

G0SPSI(36) 

sinpsi(36) 

TRIMGM 

MRA 

GH0RD(30) 

TWIST(30) 

thetzlOo) 

xa(30) 

XAG(30) 

RGMAX 

RFA 

XFA 

OPUSLD 

RIDATA 

RA(30) 

dra(30) 

MTIP 

OMEGA 

GMEAN 

FTIP(30) 

MUX 

MUY 

MUZ 

NBM 

NTM 

NBMT 

RGUST(3,3j 

CHUB(6,iD) 

RTRIGH 

XAPq(2,4,5j30) 

INGIGM 

T75 

dvbody(6) 

GONTGM 

vind(3,30,36) 

WKVIGM 

vint(3,30,36) 

interference  velocity  from  other  rotor 

WKV2GM 

gam(30,36) 

CIRG06) 

QRIGM 

save(30,36,19) 

AESIGM 

vgust(3,30,36) 

VGUSTH(3) 

gust  at  rotor  disk,  shaft  axes 
gust  at  rotor  huh,  velocity  axes 

GUSTGM 

AEROFl 


ETA(2,10,30’! 

ETAP(2,10,36) 

ETAPP(2aO,20) 

ZETA(5,30) 

ZETAP(5,30) 

DELI 


"bending  modes  at  i  =  1  to  MRA 
torsion  modes  at  r^,  i  =  1  to  MRA 


HDICH 


AEROSl 


Name :  AEROSl (ALPHA , DALPHA .COSYAW ,HACH, JPSI , IR .REVFLW , CL , CD, CM , CDR .OPTION) 
Function:  calculate  blade  section  aerodynamic  coefficients 

Corrected  Mach  number:  section  2.4.5 
Stall  model,  delayed  ex-:  section  2.4.7 
Yawed  flow,  effective  cx.:  section  2.4.6 

Calculate  2-D  airfoil  characteristics  at  effective  ©<  and  H:  section  2,4.7 
Section  characteristics  corrected  for  yawed  flow  and  stall  delay: 
sections  2.4.6,  2.4.7 

Dynamic  stall  vortex  loads:  section  2.4,7 


ALPHA 

DALPHA 

COSYAW 

MACH 

JPSI 

IR 

REVFLW 

CL 

CD 

CM 

CDR 

OPTION 


angle  of  attack  c<  (deg) 

£><c/V 
cos  -A. 

Mach  number  M 
azimuth  index  j 
radial  station  index  i 

integer  parameter:  1  if  just  crossed  reverse  flow 
boundary 

•^d 

c 

m 

^dradial 

integer  parajneter:  0  for  derivatives  of  coefficients 
in  flutter  analysis  (no  dynamic  stall  vortex  loads, 
and  calculated  data  not  saved) 


STATE(30,36,3)  AESICM 

nnr  m  a  or>  ^ 

dcdmax(30,36) 

dchmax(30,36) 

heff(30,36,3) 

aeff(30,36,3) 

DCLDS(30,36) 

DCDDS(30,36) 

DCMDS(30,36) 

MRA  RIDATA 

MCORRL(30} 

MCORRD(30) 

MC0RRH(30) 


AEROSl 


TAUL 

TAUD 

TAUM 

AEELAY 

AMAXNS 

PSIDS(3} 

ALFES(3) 

ALFRE(3) 

CLDSP 

CDDSP 

GMDSP 

OPYAW 

OPSTLL 

OPCOMP 

DEBUG 

MPSI 


RIDATA 


TMDATA 
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AEROTl 


Name:  AEROTl( ALPHA, MACH, RADIAL, OPTION, CL, CD, CM) 
Function:  interpolate  airfoil  tables 
General  reference:  section  2.4.4 


ALPHA 

MACH 

RADIAL 

OPTION 


CL 

CD 

CM 


angle  of  attack  ©<.  (deg) 
Mach  number  M 
radial  station  r/R 


integer  parameter:  if  1  calculate  c 
c^,  if  3  calculate  c^^^,  if  4  calculat 

^^2D 

”^2D 

°n»2D 


NAB 

NA(20) 

A  (20) 

NMB 

NM(20) 

M(20) 

NRB 

R(ll) 

CLT(5000) 

GDT(5000) 

GHT(5000) 


,  if  2  calculate 
e  all  three  coefficients 


AITABL 
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BODYVl 


Name:  BODYVl 

Function:  calculate  harmonics  of  airframe  motion 

General  reference:  section  5*1 *8 

DEBUG 

TMDATA 

MPSI 

MHARMF 

NBLADE 

RIDATA 

NAM 

BODYCM 

HB0DY( 16,6,10) 

RHICM 

fhub(6,36) 

AEFICM 

PHI(10,16) 

MNRICM 

KEPSI(21,36) 

TRIHCM 
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ENGNVl 


Neime :  ENGNV 1 

Function:  calculate  harmonics  of  drive  train  motion 
General  reference:  section  5*1 *9 


DEBUG  TMDATA 

MHARMF 

HPSI 

NBLADE  RIDATA 

NDH  ENGNCM 

TORQUE(36)  AEFICH 

PSID(10,6)  MNRICH 

HENG(6,10)  RHICM 

KEPSI(21,36)  TRIHCM 


MOTNFl 


Name :  MOTNFl 

Function:  calculate  rotor  generalized  forces 
General  reference:  section  5*1*7 

and  0^/^  for  trim:  section  5*2.1 


DEBUG 

HPSI 

SIGMA 

GAMMA 

MUX 

HUY 

MUZ 

chubt(i6,6) 

FHUB(6,36) 

FHUBM(6) 

QRTR(6) 

CLS 

CXS 

GTS 

CYS 

CPS 

CT 

CMX 

CMY 


TMDATA 

RIEATA 

RTRICM 


AEFICM 

QRICM 
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KOTNS 


Name •  HOTNS 

Function: 

ca] ' 

jilate  static  elastic  motion 

Genei-al  reference:  section  ^AAO 

DEBUG 

D0FA(16) 

dofd(6) 

airframe  degrees  of 
drive  train  degrees 

freedom 
of  freedom 

TMDATA 

CPRTR2 

TRIMCM 

CH'aBTl(l6,6) 

CHUBT2(16,6) 

RTRICH 

RTR2CM 

DDALF1(61 

DDALF2(6, 

1 

1 

HNHICM 

HNH2CM 

FHUBM1(6] 

FHUBH2(6, 

( 

1 

QRICM 

QR2CM 

ASPRKG(IO) 

ACNTPL(4,10) 

NAM 

BODYCM 

HENG0(2,2) 

NDM 

ENGNCH 

DELF 

DELE 

DELA 

DELR 

CONTCM 

MBl 

MB  2 

INCICM 

INC2CM 

QSSTAT(IO) 

PISTAT 

PESTAT 

HNSCM 
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BODYF 


Name :  B0DYF(LEVEL1 ,LEVEL2) 

Function: 

calculate  airframe  generalized  forces 

General  references  section  4,2.6 

LEVELl 

wake  level  for  rotor  A  and  rotor  #2:  0  for 

LEVEL2 

uniform  isiflow 

DEBUG 

HPSI 

AFLAP 

TMDATA 

GAMMA 

SIGMA 

RADIUS 

OMEGA 

0PRTR2 

reference  rotor 

TRIHCM 

VB0DY(3) 

(Jf  yf  5,) 

C»NTCM 

WB0DY(3) 

DELF 

dele 

DELA 

DELR 

DDZF 

«F 

CANTHT 

CANTVT 

BDDATA 

REULER(3,3) 

rwb(3} 

rht(3) 

rvt(3) 

velfO) 

BODYCM 

qvm(6') 

Qim6j 

QVT(6) 

SAVE(3l) 

qBDCM 

viwi(3,36) 

viHl(3,36j 

vivi(3,36) 

LMDAW1(3) 

LMDAH1(3) 

LHDAV1(3) 

WKVICM 
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BODYF 


VIW2(3,36) 

viH2(3.36) 

VIV2(3.36) 

LMDAW2(3) 

LMDAH2(3) 

LHDAV2(3) 

GWB(3)  gust  in  F  axes 

GHT{3) 

GVT(3) 


WKV2CM 


GUSTCH 


BODYA 


Name :  BO  DYA ( VWB , VHT , V VT , WMB , AFLAP , EELF , DELE , BELA , DELR , DAWB , 

FWB , MWB , FHT , FVT , ANGLES ) 

Function:  calculate  body  aerodynamic  forces 

General  reference:  section  4.2.6 

VVfB(3)  velocity  (u,  v,  w)  at  wing-body,  horizontal  tail, 

VHT(3)  and  vertical  tail;  F  axes;  ft/sec  or  m/sec 

VVT(3) 

WWB(3)  angular  velocity  (p,  q,  r) ;  rad/sec 

AFLAP  flap  angle  Sp  (deg) 

DELF  flaperon  control  (rad) 

DELE  elevator  control  (rad) 

DELA  aileron  control  ^  (rad) 

DELR  rudder  control  (rad) 

DAWB  (rad/sec) 

FWB(3)  (D/q,  Y/q,  L/q)^^;  ft^  or  m^ 

MWB(3)  (M^q,  My/q,  or  jt? 

FHT(2)  (D/q,  L/q)^Y' 

FVT(2)  (D/q,  L/q)^^;  ft^  or  m^ 

ANGLEB(6)  (cp^y^,  (3^,  ,  ^)i  deg 

CANTHT  BDDATA 

CANTVT 

LFTAW  BA DATA 

OPTINT 
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WAKECl 


Name:  WAKECI(LEVEL) 

Function:  calculate  influence  coefficients  for  nonuniform  inflow 
General  reference:  sections  3*1 *3 »  3*1 


Calculate  h  for  axisymmetric  wake:  section  3* 1*6 

Ground  effect  parameters:  sections  2.4.3,  3*1*5 

Calculate  first  blade/vortex  intersection  age  and  core  bursting 


age: 

section  3.1.7 

Wake  age  loop: 

LANDJ  =  {%-  l)*  MR*  MPSI  +  j 

JTEMJ  =  -  j 

Burst/unburst  core  radius:  section  3.1.7 

Axisymmetric  far  wake:  section  3.1.6 

Complete  C  and 

Cj^j  for  axisymmetric  geometry:  section  3.1.6 

LEVEL 

wake  analysis:  0  for  uniform  inflow,  1  for 
prescribed  wake,  2  for  free  wake  geometry 

NBLADE 

RADIUS 

ROTATE 

RRCOT 

CH0RD(30) 

MRA 

INFLOW (6) 

RIDATA 

ROTATO 

RADUSO 

other  rotor 

R2DATA 

OMEGA 

CMEAN 

RA(30) 

RTRICM 

pinter(36) 

PBURST(36} 

(rad) ;  for  rotor  (^^2 

DPSI21 

OMEGAO 

other  rotor 

RTR2CM 

BETAC 

QRIGM 

BETAS 

BETASO 

BETASO 

QR2CM 

MPSI 

DEBUG 

DEBUGV 

OPGRND 

HAGL 

debug  print  control  for  VTXL  and  VTXS 

TMDATA 
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WAKECl 


DPSI 

LSCALE 

FSCALE 

RWB(3) 

RHT(3) 

RVT(3) 

rhubO) 

RHUB0(3) 

R0FF(3) 

RSF(3,3) 

RSFO(3,3) 

KE(3) 

RFE(3,3) 

K2T 

MUTPP(3) 

KNW 

KRW 

KFW 

KEW 

RRU 

FRU 

PRU 

FNW 

DVS 

DLS 

C0RE(5) 

OPCORE(2) 

WKMODL(l3) 

0PNWS(2) 

LHW 

OPHW 

OPRTS 

VELB 

DPHIB 

D3V 

QDEBUG 

MRG 

NG(30) 

MRL 

NL(30) 


other  rotor 
other  rotor 


other  rotor 


TRIHCM 

BODYCH 

WGIGH 

WIDATA 


W2DATA 
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WAKECl 


MR 

ML 

MI 

MW 

MH 

MV 

HO 

0(3,20000) 

CNW(3,20000) 


WKClCM 
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WAKEBl 


Name:  WAKEB1(P3I,0PTICN,RBR,RBT,RB) 

F’  ":tion:  calculate  blade  position 

General  reference:  section  3»1*3 

PSI 

^  (rad) 

OPTION 

•4 

integer  parameter  controlling  calculation  of  r.  : 
if  1,  at  r_--„  and  1;  if  2,  at  circulation  stations; 
if  3*  at  inWow  stations 

RBR(3) 

rbt(3) 

^b  ’^ROOT 

r^  at  tip  (r  =  l) 

RB(3,30) 

rt  at  inflow  or  circulation  stations 
b 

MPSI 

HHARMF 

MHARM 

TMDATA 

RFA 

ZFA 

XFA 

NBLADE 

RROCT 

RIDATA 

NBM 

RA(30) 

RTRICM 

0P^KBP(3) 

MRG 

NG(30) 

MRL 

NL(30) 

WIDATA 

BETA (21, 10) 
BETAG(21) 

MNRIGM 

PSIS(IO) 

PSISO 

MMui nu 

i  '-»*  4 

ETA(2,10,30) 

ETAR(2,10) 

ETAT(2,10) 

DELI 

DEL2 

DEL3 

bending  modes  at  r^,  i  =  1  to  HRA 

bending  modes  at 

bending  modes  at  tip  (r  =  l) 

MDICM 
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VTXL 


Name:  VTXL(R1 ,R2,RP, MODEL, OPCORS, CORE, DLS , CHORD, PSI ,CPGRND,ZAGL,RTE , 

V1,V2, DEBUG) 

Function:  calculate  vortex  line  segment  induced  velocity 
General  reference:  section  3*1 •? 

Calculate:  SI  =  s^/s,  S2  -  ^2^^*  =  r^^ 

Lifting  surface  correction: 

ANGLS  =  A (deg) 

HLS  =  h  (-1.0  for  no  correction) 

RSINL  =  rsin_A_  ,  COSL  =  cos_A_,  SINL  =  sin_A. 

LLL  =  L^^,  LLS  =  L^g,  FACTl-S  = 

Image  element  in  ground  effect:  section  3*1 ‘5 


Rl(3) 

R2(3) 

RP(3) 

MODEL 

OPCGRE 

CORE 

DLS 

PSI 

CHORD 

CPGRND 

ZAGL 

RTE(3,3) 

DEBUG 

Vl(3) 

V2(3) 


r^  (at  4’) 
r^  (at  4>+aM->) 
rp  (at  P) 

integer  parameter;  1  for  stepped  vorticity  distribution, 
2  for  linear  vorticity  distribution 

integer  parameter  defining  vortex  core  type:  0  for 
distributed ,  1  for  concentrated  vorticity 

vortex  core  radius  r 

c 

for  liftiiig  surface  correction,  LT  0.  to  suppress 
:  required  for  d^^  ^  0  only 
chord  c  at  P;  required  for  d^^  0  only 

integer  parameter:  0  for  out  of  ground  effect 


7,  .  :  required  in  ground  effect  only 

AOXj 

RpE<  required  in  ground  effect  only 
integer  parameter:  debug  print  if  GE  3 


A'v  due  to  (at 
Av  due  to  Pg  (at  <^+  2i.'4^) 
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Names  VTXS(R1,R2,R3,R4,RP,H0DELT,M0  DELS , C  PCORE , CCRET , CORES , DVS , 

OPGRND, ZAGL , RTE , MDLT , MDLS , VTl , VT2 , VS 1 , VS3 , DEBUG ) 

Functions  calculate  vortex  sheet  segment  induced  velocity 
General  references  section  3*1 *8 

Image  element  in  ground  effects  section  3*1 *5 


Rl(3) 

R2(3) 

R3(3) 

R4(3) 

RP(3) 

HODELT 

MODELS 

OPCORE 

CORET 

CORES 

DVS 

OPGRND 

ZAGL 

R1E(3,3) 

DEBUG 

MDLT 

MDLS 

VT1(3) 

VT2(3) 

VS1(3) 

VS3(3) 


5 


integer  parameters  defining  trailed  and  shed  vorticity 
models  0  to  omit,  1  for  stepped  line,  2  for  linear 
line,  3  for  sheet 

integer  parameter  defining  vortex  core  types  0  for 
distributed,  1  for  concentrated  vorticity 

r  for  trailed  vorticity  (LT  0.  for  s/2) 

r^  for  shed  vorticity  (LT  0.  for  t/2) 

d^g  for  sheet  edge  test;  LT  0.  to  suppress 

integer  parameters  0  for  out  of  ground  effect 

;  required  in  ground  effect  only 

AliL 

R^;  required  in  ground  effect  only 

integer  parameters  debug  print  if  GE  3 

integer  parameters  specifying  trailed  and  shed  vorticity 
model  used 


Av^  due  to 


(at  c|>,  outside  edge) 
due  to  Fg  (at  ^  +  outside  edge) 
Av  due  to  (at  cb,  outside  edge) 

S  X 

Av  due  to  Fo  (at  cj>,  inside  edge) 
s  j 


(Av^^  =  -Avtj,  Avt^  = 

(Avs2  =  -Avsj.  Avsi^  =  -AVg^) 
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GEOHEl 


Naros :  GEOHEl  (K  ,L  ,liEVEL  iRVT  ,RWSO  ,RWSl) 
Function:  evaluate  wake  geometry 
General  reference:  section  3*1 *3 


K 

k  ( <13  =  k  ) 

L 

LEVEL 

wake  analysis:  1  for  prescribed  wake  geometry, 
free  wake  geometry 

2  for 

RWT(3) 

~T  at  tip  vortex 
w 

RWS0(3) 

at  sVieet  inside  edge 
w 

RWSI(3) 

r  at  sheet  outside  edge 
w 

TMBATA 

MPSI 

TRIMCM 

DPSI 

WIDATA 

KRWG 

KFWG 

GIDATA 

RBR(3,36/ 

BBT(3.36) 

MUTPP(3) 

DZT( 144) 

DRT(144) 

K2T 

DZSl(l44l 

DRSi(i44J 

K2SI 

DZS0(144; 

DRS0(144, 

K2S0 

1 

1 

1 

1 

WGICM 

DFWG(3,2304) 

GEOMRl 


Name:  GEOMRl (LEVEL) 

Function:  calculate  v.ike  geometry  distortion 
General  reference:  section  3*1 *3 

Prescribed  wake  geometry:  CTG  =  C^,  CTOS  =  C^*~,  TW  =  0^^  (deg) 


LEVEL  wake  analysis:  1  for  prescribed  wake  geometry,  2  for 

free  wake  geometry 


DEBUG 

MPSI 

DPS  I 

TMDATA 

TRIMCM 

NBLADE 

SIGMA 

TWIST(30) 

KHLMDA 

RROOT 

MRA 

6^^  at  i  -  1  to  MRA 

RIDATA 

LAMBDA 

LAMBDI 

interference  velocity,  due  to  other  rotor 

WKVICH 

WKV2CM 

KRWG 

OPRWG 

FWGT(2) 

FWGSI(2} 

FWGS0(2) 

KWGT(4) 

KWGSI(4) 

KWGS0(4) 

WIDATA 

CT 

CIRC(36) 

BETAC 

'DXP'PAC! 

Cp 

QRICM 

RA(30) 

MUX 

MUY 

MUZ 

RTRICM 

RBR(3,36) 

• 

WGICM 

R2S0 


-219- 


GEXDMFl 


Name ;  GEOMFl 

Function;  calculate  free  wake  geometry  distortion 
General  reference;  section  3*2 


Subprograms  required;  WGAM,  DCALC,  WCAL,  WQCAL,  VSCAL,  QSVL,  QCTL,  QVS 


DEBUG 

integer  parameter  controlling  debug 
print;  GE  1,  ixrint  D  at  ^  =  2'ir/N  each 
iteration;  GE  2,  allow  printing;  GE  3i 
controlled  by  IWGDB  and  QWGEB 

TMDATA 

MPSI 

(maximum  24,  multiple  NBLADE) 

SIGMA 

NBLADE 

RIDATA 

PHIBWG(36) 

core  burst  age  (rad) 

RTRICM 

DBV 

WIDATA 

MUTPP(3) 

DFWG(3,2304) 

WGICH 

LAMBDA 

FAGTGE 

WKVICM 

LAIfflDI 

interference  velocity,  due  to  other  rotor 

WKV2CM 

CONING 

&  (rad) 

QPICM 

circ(36) 

p/sA 

KFWG 

OPFWG 

ITEBWG 

FACWG 

WGM0DL(2) 

RTWG(2) 

nnoc'ur'/’  li  ^ 

MRVBWG 

LrawG 

ndmwg(36) 

IPMGDB(2) 

qwgdb 

DQWG(2) 

GIDATA 

D3L1 

DEB  2 

MDICM 

Name:  MINV(a,N,D,L,H) 

Function:  calculate  inverse  of  matrix 


Input: 

a(n*n) 

matrix  (destroyed) 

N 

dimension 

L(N+1) 

scratch  vector 

M(N+1) 

scratch  vector 

Output: 

a(n*n) 

A  -  inverse 

D  determinant  of  A;  0,  if  A  is  singular 


HINVC 

Name:  HINVC(a,N,D,L,H) 

Function:  calculate  inverse  of  complex  matrix 

Input: 

A(N*N) 

N 

L(N+l) 

H(N+l) 

Output: 

A(N*N)  complex  a  -  inverse 

D  complex  determinant  of  A;  0.  if  A  is  singular 


complex  matrix 
dimension 
scratch  vector 
scratch  vector 
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EIGENJ 


Marne :  EIGENJ ( N , NM , A , T , EVR , EVI , VECR , VECI , INDIC , NEI ) 

Function:  calculate  eigenvalues  and  eigenvectors  of  matrix 

Subprograms  required:  SCALEM,  HESQR,  REALVE,  CXDMPVE 
Input: 

a(n*n) 

N 

NH 
NEI 
T 

Output: 

EVR(n)  real  part  of  eigenvalues  of  A 

EVI(N)  imaginary  part  of  eigenval>.es  of  A 

VEGR(N*N)  real  paxt  of  eigenvectors  of  A 

VECI(N*N)  imaginary  part  of  eigenvectors  of  A 

INDIC(n)  if  2,  no  error;  if  1,  eigenvector  not  found;  if  0, 

neither  eigenvector  nor  eigenvalue  found 


matrix  A  (destroyed) 
order  of  matrix 

actual  first  dimension  of  arrays;  maximum  100 
0  to  calculate  only  eigenvalues 
dummy  argument  (set  to  24.  in  EIGENJ ) 


o 


DERED 

Name:  EERED(NX,NV,IX)F,C0N,A2,A1,A0,B,IX)F1  .DOFO.NAMEX.NAHEV) 

Function:  eliminate  equations  and  variables  from  system  o-^  differential 
equations 


Input: 

NX 

NV 

DOF(NX) 

CON(NV) 

A2(NX*NX} 

A1(NX*NX) 

A0(NX*NX) 

b(nx*nv) 

DOFO(NX) 

DOFl(NX) 

NAHEX(NX) 

NAt-IEVCNV) 


dimension  of  matrices 
dimension  of  matrices 

integer  vector  designating  degrees  of  freedom  to  be 
eliminated:  DOF  =  0  if  variable  not  used 

integer  vector  designating  controls  to  be  eliminated: 
CON  =  0  if  variable  not  used 

coefficient  matrices 

control  matrix 
integer  vector 
integer  vector 
vector  of  variable  names 
vector  of  control  names 


Output: 

A2  reconstructed  matrices  and  vectors 

A1 
AO 
B 

DOFO 

DOFl 

NAMEX 

KAMSV 


-224- 


QSq-RAN 


Name:  QSTRAN(MX,MX0,MX1,MV,A2,A1,A0,B0,K)F1,D0F0,NAHEX) 

Functiont  qu&sistatic  reduction  of  system  of  linear  differential 
equations 

General  references  section  6.3.2 


Input: 

A2(HX*HX} 

A1(MX*MX) 

A0(MX^»MX) 

bo(mx*mv) 

IX)Fl(MX) 

DOFO(MX) 

MX 

HXO 

MXl 

MV 

NAMEX(MX) 


coefficient  matrices 


control  matrix 

integer  vector  designating  first  order  degrees 
of  freedoms  IX)Fi(l)  =  0  for  first  order 

integer  vector  designating  quasistatic  variables: 
IX!)F0(l)  =  0  for  quasistatic 

number  of  degrees  of  freedom,  maximuji!  60 

number  of  quasistatic  degrees  of  freedom 

number  of  first  order  degrees  of  freedom 

number  of  controls,  maximum  60 

vector  of  variables  names 


Outputs 

A2  reconstructed  matrices  and  vectors 

A1 

AO 

BO 

EOFl 

NAMEX 

MX  number  of  remaining  degrees  of  freedom  (HX-MXO) 

MXt  number  of  remaining  first  order  degrees  of  freedom 


CSYSAN 


Name  s  CSYSA  N( N,MX,MX1,MV,A2,A1,A0,B0, NFREQ , FREQ , NSTEP , DOFi , FSCALE , 

NAHEX.NAMEV.NFOUT) 

Function:  analyze  system  of  constant  coefficient  linear  differential 
equations 

General  reference:  sections  7.2,  7.2.1 


N 


A2(MX*MX} 

ai(kx*mx5 

A0(MX*MX) 

B0(MX*MV) 

MX 

MXl 

MV 

IX)Fl(HX) 

FSCALE 

NAMEX(MX) 

namev(mv) 

NSTEP 

NFREQ 

FREQ( NFREQ) 
NFOUT 


calculation  control 

N  =  0  1  2  10  11  12 


eigenvalues 
eigenvectors 
check  sums 
zeros 


XXX  X  X  X 
XX  XX 

X  X 

XXX 


coefficient  matrices 


control  matrix 

number  of  degrees  of  freedom 

number  of  first  order  degrees  of  freedom 

number  of  controls 

(maximum  MX2=  2*MX-  MX1=  60;  maximum  MV  =  60) 

integer  vector  designating  first  order  degrees  of 
freedom  (zero  columns  in  A0)j  DOFl(l)  =  0  for  x^ 
first  order 

frequency  scale  f<?  cor  S~L  (in  rad/sec  to  obtain 
frequencies  in  Hz  and  times  in  sec);  there  is  no 
print  of  dimensional  eigenvalues  if  FSCALE  ^0. 

vector  of  variables  names 

vector  of  control  names 

static  response  calculated  if  NSTEP  0 

number  of  frequencies  for  which  frequency  response 
calculated;  none  if  NFREQ 0 

vector  of  frequencies  (dimensionless)  for  calculation 
of  frequency  response 

unit  number  for  printed  output 
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Output: 

LAMDA(HX2) 

MX2 


eigenvalues 

number  of  eigenvalues 

available  in  follovring  common  block: 

COMMON  /EIGVC/LAMDA(60),MX2 
CXjHjtijSX  LAMDA 


CSYSAN 
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DETRAN 


Name:  DETRAN(A,MX,HX1,MV,A2,A1,AO,BO,DOF1,NAMEX,NAME,NFOUT) 

Function:  transform  equations  to  state  variable  form 
General  reference:  section  7.1 

coefficient  matrices 

control  matrix 

number  of  degrees  of  freedom,  maximum  60 

number  of  first  order  degrees  of  freedom 

number  of  controls,  maximum  60 

integer  vector  designating  first  order  degrees  of 
freedom;  IX)Fl(l)  =  0  for  x^^  first  order 

vector  of  variable  nanies 

unit  number  for  printed  output 

Output: 

A(MX2*MX2)  coefficient  matrix 

B0(KX*MV)  control  matrix 

NAME(MX2)  vector  of  variable  names 

(MX2  =  2*MX  -  MXl) 


Input: 

A2fMX*MX} 

A1(MX*MX) 

A0(HX*MX) 

bo(mx*mv) 

MX 

MXl 

MV 

DOFl(MX) 

!'  MEX(HX) 
NFOUT 
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SINE 


Name*  SINE(W,A,ASQ,B0,MX,MX1,MV,NAHE,NAMEV,NFCUT) 
Function;  calculate  frequency  response  from  matrices 
General  reference  7.2.3 

Response  calculation:  for  last  MX  states  only 


W 

A(MX2*MX2) 

ASQ(MX2*MX2) 

bo(mx*mv) 

MX 

MXl 

MV 

NAME(HX2) 

NAHEV(MV) 

NFOUT 


frequency  (dimensionless) 
coefficient  matrix  A 

2 

coefficient  matrix  squared,  A 
control  matrix 

number  of  degrees  of  freedom 

number  of  first  order  degrees  of  freedom 

number  of  controls 

(maximum  MX2=  2^^MX  -  MXl  =  60;  maximuiri  MV  =  60) 
vector  of  variable  names 
vector  of  control  names 
unit  number  for  printed  output 


STATIC 


Name :  STATIC( A, BO, MX, MXl, MV, NAME, NAMEV , NFOUT) 
Function j  calculate  static  response  from  matrices 
General  reference!  section  7,2.2 

Response  calculations  for  last  MX  states  only 


A(MX2*MX2) 

bo(mx*mv) 

MX 

MXl 

MV 

NAME(MX2) 

namev(mv) 

NFOUT 


coefficient  matrix 
control  matrix 

number  of  degrees  of  freedom 

number  of  first  order  degrees  of  freedom 

number  of  controls 

(maximum  MX2=  2*MX  -  MXl  =  60;  maximum  MV  =  6o) 
vector  of  variable  names 
vector  of  control  names 
unit  number  for  printed  output 
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ZERO 


Namei  ZER0(A,B0,MX2,MX,MV,NX,NV) 
Functions  calculate  zeros 
General  reference:  section  7.2.4 


A(MX2*MX2)  coefficient  matrix 

B0(MX*MV)  control  matrix 

number  of  states,  maximum  60 
number  of  degrees  of  freedom 
number  of  controls 

NX  state  number  i  for  which  zeros  to  be  calculated 

NV  control  number  j  for  which  zeros  to  be  calculated 


Outputs 

LAMDAZ(MZ)  zeros  of  x./v. 

Vi  ^  J 

factor  Kj:S  x^^/v .  = 

number  of  zeros 

available  in  the  following  common  blocks 

COMMON  /EIGV2/LAMDAZ(60),K1,HZ 
COMPLEX  LAHDAZ 
REAL  K1 


-XCz.-...s.j 
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ZETRAN 


Name:  ZETRAN(Z,MZ) 

Function;  transform  matrix  for  zero  calculation 
General  reference:  section  7 *2 A 

Input: 

Z(^5Z*MZ)  matrix  A*  (A  with  x.  column  replaced  by  v.  column  of  B) 

HZ  number  of  states,  HX2 

Output: 

Z(MZ*MZ)  matrix  A.  (eigenvalues  of  which  are  the  zeros); 

the  factor  is  in  Z(MZ*MZ+l) 

MZ  number  of  zeros 

GT  0  finite  number  of  zeros  exists 
EQ  0  no  zeros,  =  Z(l) 

LT  0  X,  not  controllable  by  v. 

^  J 
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BODE 


Name  I  BODE( MX , MX 1 ,  i IV , A 2 , A 1 , A 0 ,BO , DOFl , NAMEX , NAMEV , NPLOT , NAMEXP , NAEMVP , 

NX , HV , NFO , NFl , ND,  MSCALE , NFOUT) 

Function:  calculate  and  printer- plot  transfer  function  (Bode  plot) 

General  reference:  section  7c 2.3 

A2(HX*MX)  coefficient  matrices 

Ai(MX*MXj 

A0(MX*MX) 

B0(HX*HV)  control  matrix 

MX  number  of  degrees  of  freedom 

MXl  number  of  first  order  degrees  of  freedom 

MV  number  of  controls 

(maximum  MX2=  2*MX  -  MXl  =  60;  maximum  MV=  60) 

DOFl (MX)  integer  vector  designating  first  order  degrees  of 

freedom;  DOFl(l)  =  0  for  x^  first  order 

NAMEX(MX)  vectoi  of  variable  names 

NAMEV ( MV )  vector  of  control  names 

NPLOT  frequency  response  calculation  method:  if  1,  from 

matrices;  if  2,  from  poles  and  zeros 

NA."EXP(nx)  vector  of  variable  naiaes  to  be  plotted  (inconsistent 

names  ignored) 

NAMEVP(NV)  vector  of  control  names  to  be  plotted  (inconsistent 

names  ignored) 

NX  number  of  degrees  of  freedom  to  be  plotted;  maximum  30 

NV  number  of  controls  to  be  plotted;  maximum  30 

NTO  exponent  (base  lO)  of  beginning  frequency 

NFl  exponent  (base  lO)  of  end  frequency 

ND  frequency  steps  per  decade 

(maximum  NF=  (NFl  -  NF0)*ND+  1  =  15l) 

MSCALE  magnitude  plot  scale:  if  1,  plot  relative  maximum  value; 

if  2,  plot  relative  10*^K;  if  3  plot  relative  10, 

NFOUT  unit  number  for  printed  output 
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BOEfEPP 


Name:  BOEEPP(HH,HP,NFO,NF1,ND,OPTION,NFOUT) 

Function:  printer-plot  transfer  function  magnitude  and  phase 


HM(N) 

hp(n) 

NFO 

NFl 

ND 

OPTION 

NFOUT 


transfer  function  magnitude 

transfer  function  ptese  (degrees,  -180  to  180) 

(N  =  (NFl  -  NF0)*ND+  l) 
exponent  (base  lO)  of  beginning  frequency 
exponent  (case  lO)  of  end  frequency 
frequency  steps  per  decade 

mapitude  plot  scale:  if  l,  plot  relative  maximum 
value:  if  2,  plot  relative  10**K;  if  3,  plot 
relative  10.  ^ 

unit  number  for  printed  output 
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TRACKS 


Name;  TRACKS(A2,Al ,A0,B0,MX,MX1,MV,D0F1 ,CMEGA , NAMEX, NAMSV , NPLOT, 

PER ICD, BELT, TMAX , NAMEX P, NAMEVP, NX , m' , NFCUT) 

Function;  calculate  and  printer-plot  time  history  of  time- invariant 
system  response 

General  reference;  section  7.2.5 

Calculate  eigenvalue  matrix  and  modal  matrix; 

MRED  =  H  without  unused  states  (rows) 

iiB  =  M  without  unused  controls  (columns) 


A2(HX*MX) 

ai(mx*hx) 

A0(MX*HX) 

coefficient  matrices 

30(MV*MX) 

control  matrix 

MX 

number  of  degrees  of  freedom 

MXl 

number  of  first  order  degrees  of 

freedom 

MV 

number  of  controls 

(maximum  MX2“  2'^MX  -  MXl  =  60; 

maximum  MV  =  60) 

DOFl(MX) 

integer  vector  designating  first  order  degrees  of 
freedom;  DCFl(l)  =  0  for  x.  first  order 

NAMEX (MX) 

vector  of  variable  names 

namev(hv) 

vector  of  control  names 

OMEGA 

frequency  scale  (rad/sec) 

NPLOT 

control  input  type 

1  step 

2  impulse 

3  cosine  impulse 

4  sine  doublet 

5  square  impulse 

6  square  doublet 

PERIOD 

period  T  (sec)  for  5jnpulse  or  doublet  (NPLOT  ==  3 

BELT 

time  step  (sec) 

TMAX 

maximum  time  (sec) 

(maximum  NX*NVmiAX/DELT  = 

7200) 

6) 
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TRACKS 


NAM  'XP(NX) 
NAMEVP(m') 
NX 


vector  of  variable  names  to  be  plotted  (inconsistent 
names  ignored) 

vector  of  control  names  to  be  plotted  (inconsistent 
names  ignored) 

number  of  degrees  of  freedom  to  be  plotted;  maximum  30 


NV 

NFC'JT 


number  of  controls  to  be  plotted;  maximum  30 
unit  number  for  printed  output 
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TRCKPP 


Name ;  TRCKPP( TRACS , NX , NV ,MT, BELT, NAIffiXP, NAMEVP, NFOUT) 
Functio”:  printer- plot  time  history 


TRACE(NX,NV,MT) 

NX 

NV 

MT 

BELT 

NAHEXP(frx) 

NAMEVP(NV) 

NFOUT 


array  of  time  history  traces  to  be  plotted 
nximber  of  degrees  of  freedom  to  be  plotted 
number  of  controls  to  be  plotted 
(maximum  NX*^NV  =  26) 
number  of  time  steps  to  be  plotted 
time  step  (sec) 
vector  of  variable  names 
vector  of  control  names 
unit  number  for  printed  output 
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GUSTS 


Names  GUS'^(A2,A1 , AO, B0,HX,NX1, MV, MG, DOFl.NAHEX, RADIUS, OMEGA, GRAV, 

EULER .  VSL ,  LGU3T ,  HGUST ,  NA  MEXR ,  NAMEXL ,  ML ,  NAMEXA ,  M.\  CC , 
FREQA ,RACC . NBM ,ZETA , NAIffiXB , NFCUT) 

Functions  calculate  and  print  rms  gust  response 

General  reference:  section  7.2.6 


A2(MX*MX^  coefficient  matrices 

Al(MX*MX) 

AO(MXmx) 

30(NX^‘MV)  control  matrix  (gust  in  last  MG  columns) 


MX 

MXl 

MV 

MG 


DOFl(HX) 

NAHEX(HX) 

RADIUS 

OMEGA 

GRAV 

EULER(2) 

VEL(3) 

LGUST(MG) 


mgust(mg) 

NAHEXR(3) 


number  of  degrees  of  freedom 

number  of  first  order  degrees  of  freedom 

number  of  controls  and  gusts 

number  of  gust  components 

(maximum  MX2  =  2*MX  -  MXl +  HACC  +  MG  =  6o) 

(maxim am  MG  =  3) 

integer  vector  designating  first  order  degrees 
of  freedom;  DOFl(l)  =  0  for  first  order 

vector  of  variable  names 

length  scale  R  (ft  or  m) 

frequency  scale  JT.  (rad/sec) 

2  2 

acceleration  due  to  gravity  (ft/sec  or  m/sec  ) 

trim  Euler  angles  and  (rad);  required 

for  body  axis  acceleration  only 

velocity  components  in  body  axis  frame  (divided  by 
Ah);  only  magnitude  required  (for  'x:^)  unless  body 
axis  acceleration  calculated 

real  vector  of  gust  correlation  lengths:  if  GT  0, 
dimensional  correlation  length  L  ( Xp  =  L/2V);  if 
EQ  0,  L  =  400,  used;  if  LT  0,  magnitude  is  correlation 
time  (dimensionless),  so  break  frequency  is 

U) 

real  vector  of  gust  component  relative  magnitudes 

names  of  (?ip»  r  0^^  in  state  vector  (NAMEX); 

analysis  assumes  ki&t  (^Is’  ^Is’  ^Is 

immediately  (inconsistent  names  ignored) 
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GUS-K 


HAHEXL(HL) 


ML 

NAMEXA(MACC) 


freqaCmacg) 

MACG 


names  oi  linear  degrees  of  freedom  in  state 
vector  (NAMEX)  for  dimensional  output  (ft  or  m, 
obtained  from  R);  degrees  of  freedom  not 
identified  are  angular  (degrees)  (inconsistent 
names  ignored) 

number  of  linear  degrees  of  freedom 

names  of  degrees  of  freedom  ( NAMEX )  for  which 
acceleration  calculated;  last  three  names  must 
equal  AGGB  to  calculate  body  axis  acceleration 
(all  three  or  none)  (inconsistent  names  ignored) 

accelerometer  break  frequency  (Hz),  in  same  order 
as  NAMEXA;  2/rev  used  if  FREQA  ^  0. 

number  of  accelerometers;  none  if  MAGG  $  0 


RAGG(3) 

ZETA(3,NEM) 

NEM 

NAMEXB(6+HEH) 

NFOUT 


X,  y,  z  location  of  point  at  which  bodv  axis 
acceleration  calculated  (dimensionless) 

airframe  elast  mode  shapes,  k  =  1  to  NEM;  required 
for  body  axis  acceleration  only 

number  of  airframe  elastic  modes;  none  if  NEM  ^  0; 
maximum  10 

names  of  6>pi  Xp,  yp,  Zp,  qp^^ 

in  state  vector  (NAMBX);  assumes  all  elastic  airframe 
states  are  consequtive;  required  for  body  axis 
acceleration  only  (inconsistent  names  ignored) 

unit  number  for  printed  output 
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PSYSAN 


Namet  PSYSAN(MX,MX1,A2,A1,AO,PHI,DT,NT,MT,PERIOD,DOF1,NINT,NFOUT) 

Function:  analyze  system  of  periodic  coefficient  linear  differential 
equations 

General  reference:  section  7.3 


A2(MX*MX) 

Al(MX*MXj 

ao(mx*hx) 

MX 

MXi 


DCFl(HX) 


DT 

NT 

MT 

PERIOD 

PHI 


HINT 


NFOUT 


coefficient  matrices 


number  of  degrees  of  freedom 

number  of  first  order  degrees  of  freedom 

(maximum  rDC2=  2*MX  -  MXI  =  60) 

integer  vector  designating  first  order  degrees 
of  freedom  (zero  columns  in  A0)s  DOFl(l)  =  0 
for  first  order 

time  increment;  may  vary  with  NT,  but  for  Runge-Kutta 
integration  successive  pairs  must  be  equal 

time  step  counter  (NT=0,  1,  2,  ...  MT) 

total  number  of  time  steps  in  numerical  integration; 
for  Runge-Kutta  integeration,  must  be  even 

period  T  of  the  system 

temporary  storage  of  state  transition  matrix  ^  and 
last  A;  dimension  2*MX2*MX2  for  modified  trapezoidal 
integration;  dimension  3*NX2*MX2  for  Runge-Kutta 
integration  (HX2=  2*MX  -  MXl) 

numerical  integration  method;  if  1,  modified 
trapezoidal  method,  error  order  DT**3!  if  2, 

T?n +  o  o*r'v*r\v 

unit  number  for  printed  output 


Output: 

LAMDA(MX2) 

LAMDAC(MX2) 

MX2 


roots  V  (principal  value) 
eigenvalues  \  of  ^(t) 
n’unber  of  poles 

available  in  the  following  common  block: 

COMMON  /E1GVP/LAMDA( 60) ,LAMDAC( 60) ,HX2 
COMPLEX  LAMDA,LAMDAG 
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PSYSAN 


Typical  usages 

DT  =  PERIOD/HT 
DO  1  NT  =  O.HT 
T  =  DT  *  NT 

calculate  coefficient  matrices  at  time  T 
1  CALL  PSYSAN 


DEHIAN 


Namet  IEPRAN(A,MX,MX1  ,A2,A1,A0,IX)F1,NF0UT) 

Function:  transform  equations  to  state  variable  form 
General  reference:  section  7.1 

Input: 

A2fMX*HX)  coefficient  matrices 

ai(mx*mx5 

A0(MX*HX) 

MX  number  of  degrees  of  freedom;  maximum  60 

MXl  number  of  first  order  degrees  of  freedom 

DOFI(HX)  integer  designator  of  first  order  degrees 

of  freedom;  lOFl(l)  =  0  for  first  order 

NFOUT  unit  number  for  printed  output 

Output: 

A(MX2*HX2)  coefficient  matrix  (MX2  =  2*MX-MXl) 
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HAINTB 


Name:  MAINTB 

Function:  airfoil  table  preparation 
General  reference:  section  2.4.4 

Subprograms  required:  AEROT,  AEROPP,  C81INT,  C81RD,  REDCL,  TABFIX 
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AEROT 


Name ;  AEROT( ALPHA , MACH, RADIAL, OPTION, CL, CD, CM) 
Function:  interpolate  airfoil  tables 
General  reference:  section  2.4.4 


ALJHA 

MACH 

radial 

OPTION 


CL 

CD 

CM 


angle  of  attack  ex'  (deg) 

Mach  number  M 
radial  station  r/R 

integer  parameter:  if  i  calculate  < 

®d’  ^  calculate  c  ,  if 

all  three  coefficients 


m' 


‘^*^2D 

^^2D 

^n>2D 


;  if  2 
calculate 
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AEROPP 


Names  AEROPP(CL,GD,CM,MA,AMAX) 

Functions  printer-plot  airfoil  aerodynamic  characteristics 
Calculate  ordinate  limits s 


a) 

c  =  maximum,  value  of  magnitude 

b) 

N=  [log  cl  (N  =  N- 1  if  c  <  1.) 

c) 

K  =  t_c/lO**Nl  +  1 

d) 

use  for  scale  X  =  K*  10**N 

cl(ha) 

array  of  Cjj  to  be  plotted 

gp(ma) 

array  of  c^  to  be  plotted 

ch(ma) 

array  of  c^  to  be  plotted 

MA 

number  of  angle  of  attack  values;  odd  number 

A  MAX 

maximum  angle  of  attack;  data  in  arrays  for 

to  ,  in  MA  steps 

max  max  ^ 
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.  COMPUTER  SYSTEM  SUBPROGRAMS 


The  following  computer  system  subprograms  (or  the  equivalent) 
are  required  to  determine  the  calendar  date  and  time  of  day,  which 
form  the  identification  for  jobs  and  files. 

a)  CALL  TIME(ITIHE) 

Function:  returns  time  of  day  (8  alphanumeric 
characters)  in  array  ITIME(2) 

b)  CALL  DATE(IDATE) 

Function:  returns  calendar  date  (8  alphanumeric 
characters)  in  array  IDATE(2) 

The  following  computer  system  subprograms  (or  the  equivalent)  are 
required  in  the  timing  subprogram. 

a)  CALL  SETTIM(0,,0) 

Function:  initializes  timer 

b)  ITIME  =  INTV AL( 0,0) 

Function:  returns  CPU  time,  in  milliseconds  since 
initialization 


k  CORE  REQUIREMENTS 


The  program  requires  4.04  megabytes  of  core  storage.  Of  this 
total,  1.84  megabytes  is  for  the  subprograms  and  2.20  megabytes  is 
for  the  common  blocks.  The  common  blocks  for  the  nonuniform  inflow 
influence  coefficients  (both  rotors)  require  0.96  megabytes. 
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